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High-temperature superconductors at zero doping can be considered strongly correlated two- 
dimensional Mott insulators. The understanding of the connection between the superconductor and 
the Mott insulator hits at the heart of the high-temperature superconducting mechanism. In this 
paper we investigate the zero-temperature doping-driven evolution of a superconductor towards the 
Mott insulator in a two dimensional electron model, relevant for high temperature superconductivity. 
To this purpose we use a cluster extension of dynamical mean field theory. Our results show 
that a standard (BCS) d-wave superconductor, realized at high doping, is driven into the Mott 
insulator via an intermediate state displaying non-standard physical properties. By restoring the 
translational invariance of the lattice, we give an interpretation of these findings in momentum 
space. In particular, we show that at a finite doping a strong momentum-space differentiation 
takes place: non-Fermi liquid and insulating-like (pseudogap) character rises in some regions (anti- 
nodes), while Fermi liquid quasiparticles survive in other regions (nodes) of momentum space. We 
describe the consequence of these happenings on the spectral properties, stressing in particular the 
behavior of the superconducting gap, which reveals two distinct nodal and antinodal energy scales 
as a function of doping. We propose a description of the evolution of the electronic structure while 
approaching the Mott transition and compare our results with tunneling experiments, photoemission 
and magnetotransport on cuprate materials. 



PACS numbers: 71.10.-w,71.10.Fd,74.20.-z,74.72.-h 

I. INTRODUCTION 

Since the discovery in 1986 of the high temperature (H- 
TC) superconductivity in copper oxide based materials^, 
much effort has been devoted by the scientific community 
to understand the physics of this phenomenon, but its 
key-ingredients remain still unknown. On the experimen- 
tal side, the complexity of these materials, which present 
a rich phase diagram with many competing instabilities, 
has made impossible to perform resolutive experiments. 
On the theoretical side, many theories have been pro- 
posed, but a general consensus has not been achieved yet 
because of the lack of tools capable of perform reliable 
calculations. The strongly-interacting many-body nature 
of these systems makes in fact standard techniques hardly 
applicable. The recent discovery of high-temperature 
superconductivity in a new family of materials^ (for a 
perspective see e.g. ref.— ), with a general composition 
LaFeAs0i_2;Fa;, different from the typical Cu-0 plane 
structure of all known H-TC compounds, has revitalized 
the attention on the origin of the H-TC mechanism. 

In early times, P.W. Anderson^ suggested that H-TC 
superconductivity is the result of doping a Mott insula- 
tor with a small number of carriers. The understanding 
of the way a superconductor can transform into a Mott 
insulator has been a longstanding open problem and its 
connection with a strongly-interacting superconductivity 
has been long sought. The aim of this article is to investi- 
gate the doping-driven evolution of the superconducting 
state of a doped Mott insulator in a two dimensional lat- 
tice model of strongly correlated electrons. To this pur- 
pose we employ a cluster extension of Dynamical Mean 



Field Theory (DMFT) (for a review see ref,^). 

1. Dynamical Mean Field Theory and its Gluster 
extensions to study H-TG superconductivity 

DMFT has in recent years proved very successful in re- 
vealing the physics of the Mott transition in three dimen- 
sional compounds, like e.g. Vanadium Oxide. DMFT is 
a self-consistent mean field method, which maps a lattice 
problem onto a single-impurity embedded in a bath of 
free electrons. It is exact in the limit of infinite dimen- 
sion, where the physics is purely local and DMFT virtu- 
ally provides the complete solution. Its strength stands 
in the ability of treating on the same footing high and 
low energy physics, capturing in one framework both the 
Mott insulating state as well as the metallic state, and 
hence allowing a non-perturbative investigations. A good 
summary of the physics captured by DMFT is described 
in the DMFT phase diagram. For example, single-site 
DMFT in conjunction with electronic structure methods 
has provided results in excellent agreement with experi- 
ments on numerous three dimensional compounds (with 
valence electrons in the sp shell^, 3d shell, 4d shelli, 4f 
shell^ii^ and 5f shelli^). For some recent reviews of this 
field see e.g. refjii. 

In spite of these successes in genuine strongly corre- 
lated materials, DMFT has also well understood lim- 
itations, mainly due to its purely local character. It 
does not capture for example the feedback of collective 
modes, such as magnetic fluctuations on the single par- 
ticle quantities. One consequence of this shortcoming is 
that the effective mass diverges as the Mott transition is 
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approached, while in finite dimensions exchange effects 
should remove this enhancement in the region where the 
renormalized kinetic energy is of the order of the super- 
exchange-interaction or the temperature. These effects 
are negligible whenever strong frustration, finite temper- 
ature or orbital degeneracy help justifying a local approx- 
imation. The DMFT approach, however, breaks down 
in the cuprate-based superconductors, whose physics is 
two-dimensional and which have a large super-exchange. 
As the Mott transition is approached in fact, experimen- 
tal evidence shows that the effective mass remains finite 
and that physical properties are strongly momentum de- 
pendent (see for example Angle Resolved Photoemission 
Spectroscopy ARPES result a— ) . Being a local theory, 
DMFT is not able by construction to describe a strong 
momentum dependence, hence it is not suitable to study 
real finite-dimensional H-TC superconductors. 

Recently, extensions of DMFT capable to go beyond 
the local approximation have been developed. For a re- 
view of different cluster extensions of DMFT see e.g. 
j.gf^j^ji6j2^i8_ The Cellular Dynamical Mean Field 
Theory (CDMFT— ) is a cluster method which retains 
the real space interpretation of single-site DMFT but in- 
corporates non-perturbative momentum-dependence. It 
replaces the single site of DMFT by a cluster of impu- 
rities, offering the possibility of well describing short- 
ranged spatial correlation and providing a first correc- 
tion to momentum-dependent properties. It also allows 
to naturally describe phases characterized by a spatially 
dependent order parameter, like for example antiferro- 
magnetism (AF) or d-wave superconductivity (dSC). The 
first cluster-DMFT studies of the Hubbard model in two 
dimensions were carried out on a 2 x 2 plaquette in 
ref li^iii. These groups have shown the existence of an- 
tiferromagnetism, pseudogap and superconductivity, fo- 
cusing on relatively small values of the on-site interac- 
tion (mostly U =4t), for which the Fye-Hirsch Quan- 
tum Monte Carlo (QMC)^*'-^^ method is apphcable in 
implementing the cluster DMFT. One school of thought, 
for example in refi^i^^i^, advocates the study of large 
clusters, at present possible for relatively small interac- 
tions. In this latter case, however, it is not yet known 
what aspects of the Mott physics are captured. Cluster 
DMFT implemented on the 2x2 plaquette for stronger 
values of U (the Mott regime) were the subject of several 
publications2^i2&i21i^ on the Hubbard Model. 



2. A 2x2 plaquette Cellular DMFT study 

Following the cluster DMFT studies mentioned above, 
we focus in this paper on a two-dimensional 2x2 cluster- 
plaquette. According to our viewpoint, elucidating the 
physical content of the mean field theory on a plaque- 
tte is a very important step to accomplish before pro- 
ceeding to realistic studies of the actual instabilities that 
govern the phase diagram of the model. The work in 
references22i^i^i^ on small clusters have in fact shown 



that there are several competing phases and possible 
phase-separation, which could lead to complicate pat- 
terns in real space depending on the boundary condi- 
tions or on various terms not explicitly included in the 
Hamiltonian. A proper investigation of those states in 
the framework of DMFT-based methods, requires there- 
fore an understanding of the pure phases of the simple 
cluster-plaquette. 

Implementation of CDMFT on a 2x2 plaquette for 
large values of the interaction has already revealed 
several interesting effects not present in single-site 
DMFT, indicative of a very rich physics. For exam- 
ple, publications2^i^ reveal that the approach to the 
Mott transition as a function of doping occur in a 
very anisotropic way in momentum space. In the Hub- 
bard model (relevant for the description of hole-doped 
cuprates), the spectral weight disappears much more 
rapidly in the antinodal than in the nodal region of mo- 
mentum space, where quasiparticles survive close to the 
Mott transition at a fixed temperature. Evidence for the 
formation of a pseudogap in the one electron spectra can 
be seen in other cluster-DMFT studies^i^^M. CDMFT 
allows a natural interpolation of the nodal and antin- 
odal spectral function, and it opens the way to study 
the fc-dependence of the one-electron spectra in the Mott 
insultor— . The origin of the pseudogap can be traced to 
the growth of the self energy in certain parts of fc-space, 
where eventually at zero temperature lines of poles of 
the self-energy (i.e. zeroes of the one-particle Green's 
function) appear— s^^i^. In this latter case, the Fermi 
arcs observed in ARPES are interpreted as the result of 
a binding of segments of Fermi surface and of a lines of 
zeros of the Green's function. Results along those lines 
beyond mean field theory can be found in refi^. Within 
CDMFT not only the under-doped regime is anomalous. 
At optimal doping, where the maximum of the critical 
temperature occurs, a maximum in the one electron scat- 
tering rate and the presence of power laws in the optical 
conductivity are found^. 

In this paper we study the superconducting state which 
arises upon doping a Mott insulator by implementing 
the 2x2 plaquette CDMFT with exact diagonalization^ 
(ED-CDMFT) at zero temperature (Lanczos method). 
Our work is complementary to the 2x2 plaquette 
CDMFT study of ref4i, where the continuous time quan- 
tum Monte Carlo (CTQMC) at finite temperature was 
used as impurity solver. As compared with QMC meth- 
ods, ED allows to extract in an unbiased way real fre- 
quency quantities, which can be more easily physically 
interpreted. This will turns out to be fundamental in 
interpreting our CDMFT results and making important 
connection with experiments (as we will widely present 
throughout the paper). The Lanczos method, on the 
other hand, is limited by the finite size of the system used 
to describe an effective impurity model, associated with 
the original lattice model (see the following discussion 
on the ED-CDMFT procedure). If rightly implemented, 
however, it is able to well capture the difficult physics of 
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the finite dimensional Mott transition (see for example 
references^^'^^:"''^). We will in fact show that the physics 
described by our results well compares (at least at the 
qualitative level) with the CTQMC-CDMFT results of 
reference*^, where finite size limitations are absent. 

Earher studies^'* with ED-CDMFT have compared the 
phase diagram of the two dimensional Hubbard Model 
with the one of cuprate-based materials, studying the in- 
terplay between the superconducting and antiferromag- 
netic instabilities, either in the hole and electron doped 
sides. Here we complete the work presented in a previ- 
ous short publication^, and we focus on doping-driven 
evolution of the superconductor towards the Mott insula- 
tor. In particular we show that, in a small region around 
optimal doping, novel properties appear in the super- 
conducting state which are not ascribable to the stan- 
dard Bardeen Cooper Schrieffer (BCS) theory of super- 
conductivity. The most striking phenomenon, as widely 
presented in publication^, is the rising of two distinct 
doping-dependent energy scales, which do not fit neither 
into the framework of BCS approaches nor within the 
most popular theories of H-TC superconductivity, like 
for instance the resonating valence bond theories (RVBi, 
for a recent review see e.g4^). We interpret our clus- 
ter results extracting the momentum-dependence (i.e. 
restoring the translational invariance of the lattice) of 
the electronic structure as a function of doping. Accord- 
ing to our physical picture, the Fermi liquid description 
holds at low energies in the nodal region of momentum 
space, where the d-wave superconducting gap is zero. 
We show that this fact is best described in momentum 
space introducing a periodizing scheme based on the lo- 
cal (within cluster) self-energy^^*^. On the other hand, 
in the antinodal region, where the d-wave gap is max- 
imal, besides the superconducting contribution, a fur- 
ther contribution to the one-particle gap appears also in 
the normal component at the finite critical doping. In 
this case we show that the periodization of another lo- 
cal quantity, the irreducible two-point cumulant, offers a 
more adequate description^Siil. We introduce therefore 
a mixed-periodization scheme, which was foreshadowed 
in a phenomenological approach to the transport prop- 
erties of cuprate materials in the normal state"*^ and in 
a previous CDMFT publication'*^, and confront the re- 
sulting picture of the electronic evolution as a function 
doping with photo-emission, scanning tunneling experi- 
ments and magnetotransport on cuprate materials. 



3. Set-up of the paper 

The paper consists of two main parts. The first com- 
prises sections Ill lIVI In section U after this introduction, 
we present the model and briefly explain the ED-CDMFT 
method. In section |TT] we present the raw cluster DMFT 
results, mainly stressing their evolution with doping and 
showing the appearance (at small doping) of two distinct 
energy scales. It is not however straightforward to inter- 



pret these cluster quantities in terms of physical observ- 
ables. Therefore, in section UTIl we cope with the prob- 
lem of extracting lattice quantities (which can be more 
easily compared with experimental results) from the cor- 
responding cluster ones, i.e. we present and justify dif- 
ferent periodization- methods. We show that nodal and 
antinodal regions of momentum space turn out having 
different physical properties which require different peri- 
odization schemes to be rightly described. In section HVl 
we characterize the different properties in the nodal and 
antinodal regions, according to the periodization schemes 
introduced, clearly linking them with spectroscopy exper- 
iments. 

The second part starts from section |Vl where, in or- 
der to be able to compare our results with experiments 
resolved in momentum space, we introduce a more gen- 
eral mixed-periodization scheme. This latter, interlacing 
the nodal and antinodal properties in all the momen- 
tum space, allows us to propose a description of the evo- 
lution of quasiparticle spectra in approaching the Mott 
transition. Spectra, which are directly comparable with 
photo-emission on the cuprates, are presented in detail 
in section IVll In section IVTIl we present peculiar charac- 
teristics of the spectra, the so called "kink" feature, also 
comparing with photoemission results. Finally, in sec- 
tion IVIIIl by applying a phenomenological Boltzmann 
approach on our CDMFT mixed-periodization result, we 
derive the Hall resistivity, which is a direct probe of the 
charge carries in the system. We compare its evolution as 
a function of doping with magnetotransport experiments 
and draw conclusions on a topological phase transition of 
the Fermi surface, which takes place at low doping within 
the mixed-periodization scheme introduced in section [V] 
We finally derive our conclusions in section HXl 



A. The Model 

We consider the one-band two dimensional Hubbard 
Model on a square lattice: 

, (1) 

which is universally considered a minimal description 
of cuprate-based materials^. Here Ci^a (cJa) 
struction (creation) operators for electrons of spin cr, 

"icr = c|(jCi(T is the density of electrons, /i is the chemical 
potential tuning doping and tij are the orbital hopping 
integrals. For convenience's sake we consider only the 
nearest neighbor amplitude t — 1, and a next nearest 
neighbor hopping t' = — 0.3i. We set the on-site repul- 
sion U = 12t, larger than the band-width 8t, to be in the 
Mott regime. 
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B. ED-CDMFT procedure 

Similarly to single-site DMFT^ in CDMFT the orig- 
inal model (eq. [T]) is described in terms of an effective 
action containing a Weiss dynamical field Qo (t) describ- 
ing the degrees of freedom outside the cluster (in the 
bath) as a time dependent hopping within the cluster 



•S'off 



(2) 

^ = l..Nc {Nc = 4 in the 2x2 plaquette) labels the de- 
grees of freedom inside the cluster. For the case of a 
2x2 plaquette considered in this paper, a convenient 
Nambu-spinor notation has been introduced: 



t = (J 



(Ci| , ■ • ■ , C4|, ClX, . . . , C4l) 



(3) 



With this notation the Weiss field is a 8 x 8 ma- 
trix with both normal (particle-hole) and anomalous 
(particle-particle) components'^. Physically, this action 
describes a cluster embedded in a self-consist bath of free 
electrons with dSC correlations. In the CDMFT proce- 
dure, a starting guess of the Weiss field Qo is given as 
input. Then the cluster single-particle propagator Gc is 
computed trough the effective action eq. ^ and the clus- 
ter self-energy is determined through the Dyson's equa- 
tion 



Here, 



Gc {r,T') 



Gt (r,r') 
Ft(r,T') 



F(r,r') 
-Gi {t',t 



(4) 



(5) 



is an 8 X 8 matrix, G^^^cr = {-Tc^^cr{T)cl^{{))) and = 
(— Tc^l(r)ci,|(0)), {^,i' = l...Afc label sites of the cluster) 
are the normal and anomalous cluster-Green's functions 
respectively. From the cluster self-energy Sc, we use the 
CDMFT self-consistency condition to re-compute the lo- 
cal cluster Green's function G;oc(iw„) — J2k G{K,iujn), 
where 



G{K, iUn) = iuJn + ^J■- t{K) - t.c{iU!n) 



(6) 



In eq.® t{K) is the Fourier transform of the hopping 
matrix defined on the lattice formed by the clusters and 
the sum over K is therefore performed over the Brillouin 
zone reduced by the partition in clusters of the lattice^. 
We finally re-derive a new Weiss field 
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and iterate until convergence is reached. 

In practice, as we mentioned in the introduction, in or- 
der to solve the cluster impurity problem, in this work (as 
in ref.^^'M) we use the Exact Digonalization method""'. A 



parametrized Anderson-impurity Hamiltonian describes 
the action eq. ^ and couples the cluster impurity with a 
discrete number Ni, of bath orbitals (we have fix through- 
out this work Nf, — 8, which is the limit in practice ac- 
cessible with standard computational resources) : 

Wimp = ^ E^t.^c'l^Ciy^ + ^ X! '^MT"-Pi + 

+ ^ £kaal„akcr + ^ Vkfj.aal^C^a + h.c. + 



kfia 



sup t t 



fsup 



kfj.a ^M' 



^tia^ka 



(8) 
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Here Efj^iya = -f^S^i,, S^^cr = -t5^j,±i. Under the self- 
consistency constrain eq. © and ^ , the bath-parameters 
efccr, Vk^a and Vj^^^ are determined at each CDMFT- 
iteration by fitting the Anderson-impurity Weiss field (eq. 
El) with a TVfc-pole bath function G^l'^iiujn) = iluj-E-A 



A^Y.^i^'^''~^k)-'Vk 

Vfe(2x8)=( ■'^'^L ''feMT 



Efc(2 X 2) 



-eki 



(9) 
(10) 
(11) 



The fitting is obtained via a conjugate gradient mini- 
mization algorithm, which uses a distance function^ 



/ 



eo"^'"(ic^„)-g]^7(ic^„) /u;„ 
/it/ 



(12) 



that emphasizes the lowest frequencies and it is computed 
on the imaginary frequency axis LOn = {2n— 1)tt/(3. This 
introduces an effective inverse temperature, which is a 
fitting parameter (it is not the real temperature which is 
T = in our study) and is set f3 = 300t (much higher 
than the (3 = 50< used in publication of ref.^'^) throughout 
the whole paper. This parameter determines the energy 
resolution accessed in this work (see also appendix IX)) . 
On a practical level, to start the ED-CDMFT procedure 
it is most useful to introduce a reduced parameterization 
of the bath, which enlightens the symmetries of the input- 
guess Goiii^n)- The constrain on the bath parameters can 
be then relaxed in a second step. Further details are given 
in appendix iBl 



II. CLUSTER RESULTS 

We start this section by considering raw cluster quan- 
tities, which directly output from the cluster-impurity 
solution. 
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optimal doping region. It is clear we do not intend to 
draw a quantitative parallel with cuprate-based system, 
where the typical optimal doping (unambiguously defined 
as the maximum of critical temperature Tc) is around 
Sc ^ 0.15, but we rather follow the qualitative aspects of 
the physics of these material, marking a correspondences 
with our CDMFT results on the Hubbard Model in two 
dimensions. 



B. Cluster self-energies 



FIG. 1: (Color online). Left panel: doping S — 1— < n > 
versus chemical potential /i. Right panel: cluster d-wave or- 
der parameter dOP= << c^i;C^+ii >>(t-=o) as a function of 
doping S. 



A. A d-wave Superconducting state 

As we explained above, within CDMFT a supercon- 
ducting state can be studied by allowing particle-particle 
pairing terms V^^^ak\Cqi in the impurity Hamiltonian 
eq.® (or equivalently upon a unitary rotation, pair- 
ing bath terms a^^ati, see Appendix [B|) . CDMFT-self- 
consistency condition may accept non-zero values of the 
V^™^ terms, resulting in a non-zero superconducting-pair 
Green's function F^^^{t) — Tr < c^|(r)cy| >. In draw- 
ing a parallel with the classical mean field example of 
the Ising Model-, this is equivalent to assume a non- 
zero on-site magnetization m =< Sf > as starting hy- 
pothesis in building an effective Hamiltonian, and to de- 
termine then TO self-consistently. As reported in previ- 
ous workSSj CDMFT supports indeed a d-wave super- 
conducting state in a region of doping that precedes the 
Mott transition. This is shown in Fig. [1] where we show 
the doping 5 — \— < rii > versus the chemical poten- 
tial /i and the d-wave order parameter (dOP), defined 
as ^)i,p±i(t = 0), also as a function of 5. As expected, 
5 versus /i is monotonically decreasing. The dOP has a 
dome-like shape and its sign is alternating by exchange 
of the x-y bonds on the cluster plaquette. We can lo- 
cate the maximum only with some degree of uncertainty 
around 5p ^ 0.06 < 5 < 5c ^ 0.08. The two doping 
values 5p ^ 0.06 and 5c ^ 0.08 present as two special 
points, which mark changes in the physical properties of 
the system. We will come back on these two points in 
throughout the paper. Here we stress that at 5c in par- 
ticular, our dynamical mean field result seems to branch 
two different lines of solution. This is evident either in 
the 5 — fi and the dOP plots. In mean field approaches 
this behavior may be the signature of a phase transition. 
It is therefore intriguing that we find such a behavior 
close to the maximum of the dOP dome. 

Following the nomenclature typical for cuprate materi- 
als, we will call hereafter the region 5 < 5p under-doped, 
5 > 5c over-doped and 5p ^ 0.06 < 5 < 5c ^ 0.08 



It is worth to investigate in detail the cluster-outputs. 
The typical output of the CDMFT-scheme is a cluster- 
self-energy (eq. IH), which can be expressed as a 2x 
Nc X Nc {Nc = 4 for the 2x2 plaquette) matrix with nor- 
mal and anomalous components in the Nambu notation 
introduced in eq. ([3]). We first look at the normal com- 
ponents, which arc shown on the Matsubara- frequency 
axis in Fig. [2l The real parts are displayed on the left 
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FIG. 2: (Color online). Normal components of the cluster 
self-energy vs the Matsubara frequency ujn as a function 
of doping 5. In the left column the real parts are displayed, 
in the right column the imaginary parts. From the top row 
to the bottom we show the local En, the next neighbor E12, 
and nearest next neighbor S13 self-energies. 



and the imaginary parts on the right column from high 
(in our solution) doping 5 — 0.14 (the over-doped side) 
until close to the Mott transition for 5 — 0.02 (in the 
under-doped side). From the top row to the bottom we 
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FIG. 3: (Color online). Eigenvalues of the normal component 
of the cluster self-energy E^,y vs the Matsubara frequency ujn 
as a function of doping S. In the left column the real parts are 
displayed, in the right column the imaginary parts. From the 
top row to the bottom we show the self-energy corresponding 
to the points k = (0, 0), (0, n) and (tt, tt) in the first quadrant 
of the Brillouin Zone (see section IIII|I . 



have the local self-energy Sn, the next neighbor S12 and 
the nearest next neighbor E13. At high doping, the local 
components is dominant. We expect in this region that 
the single site DMFT is already a good approximation. 
By decreasing S, however, also the non-local components 
S12 and S13 grow considerably, becoming comparable 
with the local one. This effect is totally missed by single 
site DMFT and it is only captured by using a cluster ex- 
tension. The growing of the non local components of E^i^ 
determines the physical properties of the system in ap- 
proaching the Mott transition in a fundamental way, as 
we explain in the following. The imaginary parts (right 
column) display a low-energy Fermi-liquid-like behaviour 
ImE ^ aujn ■ The slope a is connected with the quasipar- 
ticle residuum (see the following). Here we notice that, 
with respect to the other components, the energy range 
of linear behavior in the nearest next neighbor lmSi3 
is narrower and that its slope a also changes behavior 
in passing through the doping Sc ^ 0.08 (it grows for 
S < 6c, decreases until becoming negative for S > Sc), 
while it also changes sign at Sc ~ 0.06. This behavior 
is the only mark we find in the cluster self-energy of the 
special nature of the points 6c and 6p that we identified 



in the Fig. [TJ The local and next neighbor components 
show instead a smooth continuous behavior as a function 
of doping i5. 

As in ref.^'''^^, it is convenient to look at the eigen- 
values of the cluster-self-energy matrix, which can be di- 
rectly related to the corner points of the quadrant of the 
Brillouin Zone (as we will explain in more detail in the 
following section). A convenient way is to re-express the 
Nambu-spinor notation eq. ([3]) by grouping two by two the 
up and down spin construction/destruction operators on 
each cluster-site: 



= (4|,Ci;) . . . , (4^,C4x) 

The cluster self-energy matrix assumes the form: 



(13) 



S2 Sij, 



Sly \ 

S2 

So Sia 
S12: So / 



(14) 



where the 2x2 matrices are function of cluster elements 
y 



So 
S2 



^ — 1^1 ,CT ■ 




Sot 
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S2t 
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2i 



S12; 





^ano \ 




-SiiJ 








-Sii 



(15) 



We follow the procedure used to diagonalize the cluster- 
self-energy matrix in refi^: 



/Soo Q \ 



Eott 

V S^o/ 



(16) 



The diagonal elements are linear combination of the orig- 
inal cluster self-energies matrices: 



Soo = So + Si + S2 

Sott = So + Sq — S2 

SttO = So — Sa — S2 

Stttt — So — Si — S2 



For convenience's sake we have defined 



Si 



Si| 

-Su 



-'ano 





(17) 



(18) 



Notice that from eq. (fTT)) the anomalous self-energy ap- 
pears only in the (0, tt) and (tt, 0) components, where 
the dSC gap is expected to open. The 2x2 eigenvalue- 
matrices of the cluster-self-energy are interpreted as de- 
scribing the four momentum-space points (0,0), (0, tt), 
(tt, 0) and (7r,7r) (hence the choice of the labels), as we 
will explain in detail in the following section Hill 

The normal components of the eigenvalues of are 
shown in Fig. [31 in a fashion similar to Fig. |21 The 
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FIG. 4: (Color online). Cluster quasiparticle residua, associ- 
ated with the eigenvalues of the cluster self-energy matrix, as 
a function of doping S. 



first striking difference with respect to the Fig. [2] is that 
the eigenvalue Sqo changes little as a function of dop- 
ing (either the real and the imaginary parts). Sqtt and 
E^^ both grow in reducing doping S, but the growth in 
Etttt is on order of magnitude bigger. This eigenvalue 
is indeed the one that drives the system into the Mott 
insulating state as (5 ^ 0, as evident from the big val- 
ues of both the real part and of the slope of the imagi- 
nary part at small doping. This behavior are very sim- 
ilar (i.e. it appears as a smooth continuation) to the 
one displayed by the self-energy of the normal state ED- 
CDMFT study published in ref,^. With respect to this 
latter study however, here we display the normal com- 
ponents of a superconducting solution (while in refi^ no 
superconductivity is allowed) and the energy resolution 
achieved {(3t = 300) is an order of magnitude smaller 
than in ref.^^ {(3t = 32). This allows us to extract the 
very low-energy properties. In particular, as expected by 
the fact that E has to be negative in order to respect ca- 
suality, we can observe in Fig. [3] that the slope a of the 
ImE ~ aujn is always negative. The fact then the w„ — s- 
behavior of ImE is linear indicates that these eigenval- 
ues components have Fermi liquid properties. Following 
a standard Fermi liquid approach, it is instructive to de- 
fine cluster quasiparticle residua Zx = (1 — dT,x / dujn) 
with X = (0, 0), (0, tt), (tt, tt) , even if these quantities 
have a real physical meaning only in correspondence of a 
real Fermi surface (i.e. Zx has to be interpreted here as 
renormalized cluster quantity). In the standard picture of 
the Mott transition^, the quasiparticle residuum Z ^ 
as doping is reduced S ^ 0. From Fig. |4] we see that 
only Z(7r,,r) appears to display this behavior, while Zjq.o) 
and -^(o.tt) clearly extrapolate to a non-zero value. A last 
remark concerns once again the special critical doping 
dc, where remarkably the -2^(0, tt) shows a clear change in 
behavior. For 6 > 6c it closely follows Z(^_„y decreasing 
with doping as expected in the standard Mott transi- 
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FIG. 5: (Color online). Left panel: real part of the anomalous 
component of the cluster-self-energy Eano (the imaginary part 
is negligible) vs Matsubara frequency ujn for different doping 
S. Right panel: ReEano(uJn — » 0) versus doping 5. 



tion picture. At Sc, however, Z(o.7r) departs from Z(7r.7r), 
and shows a behavior more similarly to ■2(o,o)- This once 
again indicates a quick change in the physical properties 
of the system in correspondence of 6c- We will come back 
later to discussing the physical interpretation of these ob- 
servations on the cluster quasiparticle residua (see section 

m. 



We now turn to the anomalous component of Ec, which 
we show in the left side of Fig. [5] on the Matsubara 
frequency w„. As resulting output of our CDMFT so- 
lution, only the real part of the nearest-neighbor com- 
ponent is appreciably non-zero on the Matsubara axis, 
and it assumes a d-wave sign bond-alternating value on 
the cluster-plaquette. ReE^no starts assuming apprecia- 
bly non-zero values in the over-doped side i5 = 0.14, and 
grows for decreasing doping, until 6c ^ 0.08 is reached. 
In the under-doped side however at (5 < ~ 0.06 the 
curves change tendency, and they decrease by further de- 
creasing doping. The qualitative change in the behavior 
is better enlightened by looking at the w — *■ limit: in 
the under-doped side the curves reach uj — with a fi- 
nite slope, while in the under-doped side (after 6 < 0.08) 
the slope has a smaller value. In right side of Fig. [5] we 
show the Wn ^ extrapolated value of ReEa„o as a func- 
tion of doping 6. In first approximation this value can be 
related to the superconducting gap (as we will explain 
in detail in the section IIIIl see eq. [50]) . The behavior 
of the ReEano as a function of doping shows therefore 
to be non-monotonic and roughly tracks the behavior 
of the order parameter (see Fig. [1]), similarly to stan- 
dard BCS theory. This result is fundamental and it is 
in striking contrast with slave-boson resonating valence 
bond theories^i^i^, where the amplitude associated to 
a particle-particle paring channel is a monotonically de- 
creasing function with doping and it has its maximum 
close to the insulating transition. In section IIVI we will 
better discuss the physical consequences of these results, 
making connection with recent spectroscopy experiments 
on cuprate-based materials. 
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With the ED-CDMFT procedure it is straightforward 
to analytically continue on the real axis^ (differently 
from other computational methods like e.g. Quantum 
Monte Carlo, which require further approximate meth- 
ods). The cluster Green's function Gc(wn) and the clus- 
ter self-energy ScCwn) (via the Dyson's equation [4]) , are 
expressed in a pole expansion form containing terms of 
the form hj {uj„—a). To go from the Matsubara to the real 
axis it is then enough to substitute iw„ ^ uj + irj, where 
77 is a small broadening introduced to display poles (for 
more details see appendix |^ . The price to pay comes 
from the discreteness introduced by truncating the bath 
in the impurity model (eq. [8]) with a finite number Ni, 
of orbitals. This reproduces continuous functions, like 
Gc(tf„) and Sc(t^n), through a finite number of poles. 

In Fig. [S] we display on the real axes the local den- 
sity of states N{u!) = - [Gc{i^)]ii, obtained from the 
cluster- impurity output (see eq. |4]), for different doping 
6 (from top to bottom), using a small imaginary broad- 
ening i77 = 7i X 10~^. As we will explain in detail in 
the following section, in a d-wave superconducting state 
a linear in uj "V" shaped density of states is expected for 
w — > 0. This behavior cannot be captured by the dis- 
creteness of our ED-CDMFT solution for small ui, hence 
we have rather a "U" shape for w — > given by the 
broadening rj. The shape of N{uj) and the energy scales 
of the superconducting gap can however be estimated by 
the location and intensities of the peaks. So, in looking 
at N{uj) from the over-doped side (top row, left panel 
for 5 = 0.16) to the under-doped side (bottom row, left 
panel for 5 — 0.02) of the phase diagram, we can make 
interesting observations: 

1. N{lu) is asymmetric function around a; = at low 
doping ((5 < (5c)- A rather symmetric shape is in- 
stead observed around optimal doping 6c ^ 0.08 > 
^ > (5 ~ 0.06 (in agreement with previous cluster 
DMFT resuhs of ref.^i) . 

2. the total superconducting gap, which can be eval- 
uated by measuring the distance of the spectral 
peaks from a; = 0, is increasing by reducing doping 
6. 

To better elucidate these observations, we have measured 
the distance from the Fermi level of the left (red circle) 
and right peaks (green square), and displayed them as 
a function of doping S (bottom row, right panel). In 
the over-doped side {6 > 5c) the peaks are equally dis- 
tant from UJ — 0, i.e. the superconducting d-wave gap is 
symmetric at low frequency, as expected in a standard 
d-wave BCS theory. At optimal doping Sp < S < Sc, 
however, the asymmetry in the gap strikes in, and it is 
present in all the under-doped side. This is in agreement 
with the cluster DMFT results obtained in ref4i and 
it is in nicely agreement with Scanning Tunneling ex- 
periments (STM--) on cuprate materials. The fact that 
N(uj) is most symmetric at optimal doping 5p < 5 < 6c 
(as evident from the red circles which display a change 
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FIG. 6: (Color online). Local density of state —^Gii{lo) 
obtained from the cluster-impurity solution (eq. [8]) for var- 
ious doping 5 (from top to bottom). The broadening used 
for display is 77 = 7 x 10~^. A measure of the energy scale 
associated with the one-particle gap (always present in the 
superconducting state) is given by the left (red circle) and 
right (green square) distances of the ED poles from the Fermi 
level UJ = 0. These are displayed as a function of doping 5 in 
the right bottom panel. 



in the slope as a function of S) is in line with the obser- 
vation carried in Fig. [T] on the special nature of this re- 
gion (consistently e.g with the avoided-quantum-critical- 
point scenario proposed in the CDMFT study of refi^). 
The black dashed line is the average of the right and left 
gaps, and the fact that it is increasing by reducing doping 
(similarly to the predictions of resonating valence bond 
slave boson theories^'"*^'^^) shows the important result on 
the presence in our solution of another energy scale, dif- 
ferent from the energy scale marked by the behavior of 
KeY,anoi^n ^ 0) (which is instead decreasing at low dop- 
ing {S < Sp), see right panel Fig. [S]). The interpretation 
of these two energy scales brings to interesting physical 
insights, which can be put in relation with recent exper- 
iments on the cuprate materials^'^'^'^'^'^'^'^. This 
has been presented in a previous short publication^ and 
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it will be re-proposed more in detail in the following sec- 
tions. 



III. MOMENTUM DEPENDENT QUANTITIES: 
PERIODIZATION PROCEDURES 



We want now to interpret the cluster results we pre- 
sented in the previous section in terms of physical observ- 
ables, which could be possible related to experiments. 
The relevant information is typically embodied in the 
one-particle Green's function, which in a superconduct- 
ing state can be conveniently written in a Nambu-matrix 
notation: 



^a7io 



(fc,a;) uj + Ck 



(19) 



Here — tk — fJ- is the free band dispersion of our model 
(see equation [1}. In order to determine the Green's func- 
tion we need therefore to determine the momentum de- 
pendent self-energy from the cluster solution, i.e. we 
need a periodization scheme. In previous wor k^^i^^i^'^ 
various periodization schemes have been proposed. The 
idea consists in determining the most local quantity 
Wr, which can be captured within the dimension of 
the cluster-impurity, and construct its truncated Fourier 



expansion 



15. 



(20) 



Smaller the neglected Fourier coefficients are (for \fi — 
v] > ViVc), compared to the cluster W^ct(|m ~ niore 
the fc-dependent quantity Wo-(fc) is well approximated. In 
the following we present two possible cluster-quantities 
that can be adopted to construct the fc-dependent self- 
energy, showing in which cases they can be considered 
good local quantities. 



A. Self-energy S-periodization 

The cluster self-energy (eq. [H)) is a natural candi- 
date. It is convenient, for the discussions in the following 
sections, to recast formula ((20)) in terms of the cluster 
eigenvalues (eq. [T| 



(21) 



X 



where ^^(fc) is a 2 x 2 matrix containing normal and 
anomalous components: 



S^f(a;) Sa„o(w) 



(22) 



We notice that with this formula the cluster eigenvalues 
Sx are directly related to the corner points in first quad- 
rant of the Brillouin Zone X =(0,0), (0,7r), (7r,0) and 
(tt, tt) (and we justify the notation introduced in eq. I16p . 
jx{k) are positive functions, such J2x lx{k) = IVfc^i^: 



700 = 




+ COS kx 


+ cos ky 


+ COS kx COS ky) 






— cos kx 


— cos ky 


+ COS kx COS ky) 


707r = 




+ cos kx 


— COS ky 


— COS kx COS ky) 


7x0 = 




— cos kx 


+ COS ky 


— COS kx COS ky) 



(23) 

We remark that, by construction, with this procedure we 
assume that the system is a simple Fermi liquid. We have 
in fact shown in Fig. [2]that the normal components of the 
eigenvalues of the cluster self-energy have Fermi-liquid 
behavior (ImEjf — > for w ^ 0). Their simple linear 
combination extends this property to all the fc-space. The 
anomalous component of the lattice self-energy Tiano{k) 
turns out to have a d-wave shape 



^anoik) — 2X](jno [coskx COS ky) 



(24) 



in agreement with the symmetry of the superconductive 
gap measured in experiments on cuprate s^^' ^"^ 



B. Cumulant Al-periodization 

In a normal state study of the two-dimensional Hub- 
bard Model of refi^i^, it has been shown that a more 
suitable local quantity to describe the Mott transition is 
the two point irreducible cumulant M , which arises from 
the atomic limit by perturbatively expanding the hop- 
ping term t in Hamiltonian ([T]). It is simply related to 
the normal-component lattice self-energy E^°'' by: 



(w) 



(25) 



In the cluster-impurity we have 2 Nc x Nc cumulant re- 
lations (2 is for the spin degeneracy), conveniently rep- 
resented by a Nc x Nc cumulant-cluster-matrix Mac- 



MTJH = i^Lo + ^l)l 



Snor 



(26) 



where 1 is the Nc x Nc identity matrix. The eigenvalue 
of the cumulant matrices are straightforwardly related to 
the eigenvalues of the self-energy matrices (eq. [T6|) : 



{iLu + fi - Y.'^""- ) 



(27) 



where the notation is again X =(0,0), (0,7r), (7r,0) and 
(7r,7r). As in the eq. (HU, we obtain the lattice cumu- 
lant A4^°^ by periodizing the eigenvalues of the cluster 
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FIG. 7: (Color online). We compare the nodal and antinodal 
quasiparticle peaks obtained from the normal component of 
the superconductive ED-CDMFT-solution at zero tempera- 
ture and the QMC-CDMFT normal state solution at much 
higher temperature (l/T = 75). The functions on the real 
axis with QMC are obtained with the maximum entropy 
method. To display the ED-CDMFT curves the slightly uj- 
dependent broadening r\(u)) = 0.075 -I- uP' was used (see dis- 
cussion in appendix . 

cumulants: 

Mr ^Y.^Tlx{k) (28) 

X 

Inverting eq. (|25p . we finally obtained the normal- 
component lattice self-energy J]™'" in the M.- 
periodization. 

C. Nodal and Antinodal dichotomy: S-vs 
TVf-periodization 

It is not trivial to decide whether it is better to peri- 
odize the self-energy S or the cumulant M.. The choice 
could be strongly dependent on the physical properties of 
the system, which are not a 'priori well known. And it is 
not trivial either understanding a priori to which extend 
the two approximated schemes could be able to describe 
such properties. In this case, we rely on experimental 
results on cuprate materials, like e.g. the already men- 
tioned ARPESi^!^, to fix a physically reasonable starting 
hypothesis. It is a well established experimental fact that 
approaching the Mott insulator the normal state Fermi 
surface (measured at temperatures above Tc) breaks 
up, displaying well defined quasiparticles in the regions 
(nodes) close to the centers of the quadrants of the Bril- 
louin Zone fc ~ (±^,±^), while quasiparticle disappear 
in the strong scattering regions (anti-nodes) close to the 
corners of the quadrants k ^ (0,±7r) and (±7r,0). Be- 
sides photo-emission experiments, measures of transport 
properties show anomalous (non Fermi liquid) power-low 



exponents in the temperature-dependence, especially in 
the under-doped regime of cuprate materials. In this 
case, a series of phenomenological approaches, which try 
to explain these experimental observations in the frame- 
work of Boltzmann theory^i^iSi^, have been based on 
dividing the momentum space in regions of high quasi- 
particle scattering rate, hot spots around the antinodal 
points, and regions of low quasiparticle scattering rate, 
cold spots around the nodal point, where quasiparticles 
have a much longer lifetime. The simple idea underlying 
this choice is that the nodal region behaves as a standard 
Fermi liquid, while the system in the antinodal region is 
like an insulating state. 

This nodal/ antinodal dichotomy fits into the CDMFT 
frame as the natural path taken by the system to ap- 
proach the doping-driven Mott transition2&i21i2&dl. This 
is a general property independent of the periodizing 
scheme adopted. Some important differences however 
hold between the S and the cumulant M. periodizations, 
which we will illustrate in the following. 

In Fig. [7] we show for example the spectral func- 
tions A{k,u}) = -ilmG^°''(w) in the nodal fc ~ (f , f ) 
(black continuous line) and antinodal k ~ (0, tt) (red 
dashed line) points of momentum space, obtained via M- 
periodization for the case U /t = 12, t' = — 0.3t and 5% 
doping. We confront our zero temperature ED-CDMFT 
result with a QMC-CDMET^i^s at much higher tem- 
perature (r = l/75t), where the system is in the normal 
state. In the ED case we have extracted the normal part 
of the superconducting solution (setting S^no = in eq. 
[TO]) , interpreting it as a low-temperature normal state 
parent of the high-temperature QMC-solution. This is 
far from being a trivial statement, as the normal com- 
ponent of a superconducting solution is not generally a 
normal state solution. Our aim in this picture is how- 
ever to present a qualitative comparison between two 
very different impurity-solving methods in two very dif- 
ferent regimes to show the generality and solidity of the 
nodal/antinodal dichotomy concept. This figure serves 
also to compare results of maxent analytic continuation 
of Fye-Hirsch QMC scheme with the ED results at low 
temperatures. The qualitative agreement is reasonably 
good, supporting the observation we made above on ex- 
perimental results and the dichotomy nature of the phys- 
ical properties in this system. 

We observe in particular in Fig. [7] that the antinodal 
quasiparticle peak is very broad in the high-temperature 
QMC solution, denoting a short hfetime of quasiparticles, 
while it is sharper in the low-temperature ED solution, 
much more than the narrowing due to the different tem- 
perature (we have scaled the heights of the peaks taking 
the high of the nodal quasiparticle peak as reference to 
fix the scale between the ED and QMC curves). This 
goes in the direction of ARPES experimental observa- 
tions, which show a sharpening of the antinodal quasi- 
particle in going from the normal to the superconductive 
state^^ by decreasing temperature. We finally point out 
that with the A^-periodization used in Fig. [3 in the 
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antinodal point the quasiparticle peak shifts to negative 
energies, opening a pseudogap. This is also in line with 
the above mentioned ARPES observations. On the con- 
trary, by using the S-periodization the antinodal peak re- 
mains always (for all dopings) at the Fermi level, but with 
reduced spectral weight compared to the nodal poinfe2&. 

All these observations, either from experimental facts 
and from the output of our method, point towards a di- 
chotomy of the nodal and antinodal regions of momen- 
tum space, which show different coherence energies. The 
nodal region has the highest coherence scale and sharp 
quasiparticles. In the antinodal region a pseudogap opens 
in the spectrum, and quasiparticles, if present at the gap- 
edge, are more broad and incoherent. It is therefore natu- 
ral to assume that the Fermi-liquid nodal region is better 
portrayed by the S-periodization (which as we mentioned 
above describes a Fermi liquid by construction), while 
the antinodal insulating region is better described by the 
Al-periodization (as shown in Fig. [7]). 

Taking this last assumption as starting point, we com- 
pare the results of the S and A^-periodization, stress- 
ing virtues and defaults in the frame of the physical ob- 
servations carried over above. A stringent test is given 
by reconstructing the local density of states N{u!) — 
— i ImG'^(w)ii (where Gf.{uj) is obtained by eq. [T9)) . 
which can be compared with the local Green's function 



Gi 



{Gc)ii, obtained directly in the cluster- impurity 



solution. This test is presented in Fig. [H First, it is once 
again instructive to separately study the normal compo- 
nent of the superconducting state by setting Y,ano = 
in eq. ([TO]) . We focus the attention on a small-doping 
case 5 — 0.05 close to the Mott transition point. The 
top panel displays a full energy range which includes the 
lower and upper Hubbard bands. The cluster-impurity 
density of states N{uj) is represented by the continuous 
black line, the Al-periodization result by a red dot-dash 
line and the E-periodization result by a green dash line. 
Already at small doping, in the region labeled (A), the 
formation of a Mott gap is visible in the cluster-solution. 
The Al-periodization is capable to describe this part of 
the spectrum quite well (matching the continuous black 
line of the impurity result). On the contrary, the E- 
periodization creates artificial states in the Mott gap. 
This is true also at low energy (region labeled B). A 
close up on the Fermi level {—t < uo < t) is displayed 
in the middle panel (with again Y^ano — 0). The M- 
periodization reproduces the formation of a low-energy 
pseudogap in the normal component of the spectra, as 
already remarked in previous normal-state studies^i^, 
while, once again, the E-periodization introduces artifi- 
cial states that fill the pseudogap. This simple test there- 
fore evidences the failure of the E-periodization in well 
describing the high and low energy normal-component 
of Green's function, in the regions of momentum space 
where a pseudogap is present, like in the antinodal k- 
point of Fig. [71 The scenario is different if we look at 
the low-energy superconducting density of states N(lij) by 
restoring Eq„o 7^ in eq. p^ . as presented in the bot- 
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Top panel: reconstructed lo- 
G^'^iuo) for the normal com- 
ponent (Eano = 0) of the Green's function at small doping 
((5 = 0.05). G^°^ (uj) has been obtained by either the cumulant 
A^-periodization (red point-dash line labeled M-scheme) and 
by the E-periodization (green dash line labeled E-scheme), 
and confronted with the local density of states of the cluster- 
impurity solution (black continuous line). Middle panel: 
close-up at the Fermi level. Bottom panel: the local den- 
sity of states of full superconducting solution (Ea„o 7^ 0) se- 
lects in the fc-summation only the nodal points. A low-energy 
confront with the cluster-impurity solution shows a different 
matching of the cumulant and self-energy periodizations with 
the cluster-impurity result. The display on the real axis of 
the Green's function has been obtained by introducing a uo- 
dependent broadening irjiuj) (see appendix \K\ for details). 
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torn panel of Fig. [8l Here it is the S-periodization line 
(green dash line) which better portrays the cluster spec- 
trum (black continuous line). The A^-periodization this 
time introduces spurious states close to the Fermi level 
{u! = 0). This is not contradictory with respects to the 
result presented above on the normal component of the 
system (top and middle panels). By adding the d-wave 
superconducting gap in fact, we selects at low energies 
{uj 0) only the region of momentum space close to the 
nodal points, where the gap is zero. These points domi- 
nate in fact the sum — ^ J^k IinGj,((jj)ii. N{uj) for a; ^ 
is therefore a direct probe of the nodal point behavior 
only (while without superconducting gap, i.e. Sano = 0, 
all fc-points eventually contribute in the summation at 
low energies). 

We clarify this last statement. Let's assume as starting 
point that at the nodes quasiparticles are well defined on 
the full range of doping. We can extract from eq. 
the low energy (a;„ — *■ 0) Green's function: 

ReEfe(c<j„) - Rel]fe(0) 
ImSkiuJn) ~ (1-Zj:^)w„ 



and 



+ (28) 



lUJn + Ek 



For convenience's sake, we have enlighten the quasipar- 
ticle dispersion 



6 = Zk [a+ReEfe(0)] 



Ek = ^J\U^ + \Zk^ano{k,OW 



(30) 



The nodal point is the only one gapless at low energies 
and the quasiparticle spectrum can be linearized (eq. I30p : 



Eh 



\Zk„ 



o{knod)\'^ 



ll2 /r-2 



(31) 



Viiod — |Vfc^fc| is the quasiparticle Fermi velocity perpen- 
dicular to the Fermi surface, Vano = Z„od\'^k^ano{k)\ 
is parallel to the Fermi surface. After analytic contin- 
uation ujn ^ Lo + iS, it is now easy to calculate the 
low energy behavior of the one-particle density N{uj) — 



iV(w) 



1 Zjiodj 



(32) 



i.e. Niuj) is linear in frequency close to the Fermi energy 
and the slope is uniquely determined by the quasipar- 
ticle nodal velocity Vnodi by the nodal derivative of the 
superconducting gap Vano and by the nodal quasiparticle 
residuum Zk^^^. In other words, the low-energy density 
of states in the d-wave superconductor is a direct mea- 
sure of the spectra at the nodes. This result explains 



the low-energy spectrum of the cluster density of states 
(black-continuous line in Fig. [5]), which is roughly linear 
for ^ 0, in the limit of the energy-resolution given by 
ED impurity solver. 

The fact that here the E-periodization better portrays 
the cluster result is consistent with the hypothesis of 
Fermi-liquid behaviour at the nodal points of momen- 
tum space, well described by the self-energy. The results 
of the normal component of the system suggest instead 
that in other regions of momentum space (at least close 
to the Mott transition point), a pseudogap opens in the 
spectrum and quasiparticles die in the anti-nodes. In 
the latter regions a cumulant A^-periodization is more 
appropriate than the S-periodization. 



IV. PROPERTIES OF THE NODAL AND 
ANTINODAL POINTS 

In this section we follow the observations presented 
above on the nodal/antinodal dichotomy and study the 
physical properties of the system in the nodal and antin- 
odal points, completing in detail the work presented in 
ref4^. In order to be able to make contact with the ex- 
periments, we need to extract the fully momentum de- 
pendent Green's function G'fe(cj) (see equation [T9|). which 
can be generally related to the response functions. We 
employ therefore a periodizing procedure to extract the 
/c-dependent normal and anomalous components Sfe(cL;) 
and Y,ano{k,oj) of the fc-dependent self-energy. As dis- 
cussed above, it is reasonable to start by assuming a d- 
wave shape of the superconducting gap, which is natu- 
rally obtained in our scheme by periodizing the anoma- 
lous component of the cluster self-energy via eq. (|24p . 
At the nodes, we periodize also the normal component of 
the self-energy via eq. (|2ip . This guarantees in particular 
Fermi liquid properties in the nodal region of momentum 
space. At the anti-nodes instead the insulating properties 
of the normal component of the system (like the forma- 
tion of a pseudogap at low energy and of the Mott gap 
at higher energies, see Fig. [5]) are better portrayed by 
periodizing the cumulant according to eq. (|^ . 

By adopting this procedure we can first obtain quasi- 
particle spectra in the nodal and antinodal points of mo- 
mentum space, as it is shown in Fig. [51 In order to 
define the k vectors in the nodes and in the anti-nodes, 
we follow an operative definition, similar to the one used 
in the ARPES experiments of refi^. The nodal point 
knod [antinodal point fcanorf] is chosen as the one having 
the sharpest quasiparticle peak in moving on the path 
(0, 0) (tt, tt) [(0, tt) — > (tt, tt)] of momentum space. The 
vectors knod and kanod as a function of doping 5 are shown 
in panel A of Fig. [51 We notice that kanod is a monotonic 
decreasing function of the 5. This is intuitively expected 
in a standard Fermi liquid, where the approach to the 
Mott transition at (5 = is accompanied by an increase 
of the volume enclosed by the Fermi surface in momen- 
tum space, which is proportional to the density n = 1 — 5 
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FIG. 9: (Color online). Panel A: nodal fc„od and antinodal 
kanod positions of the quasiparticle peaks in the first quadrant 
of the Brillouin Zone as a function of doping 5. Panel B: 
Nodal peaks are found always at the Fermi level and decrease 
with doping 5 in approaching the insulator. Panel C: At 
high doping {5 > 5c ^ 0.08) antinodal peaks are also found 
at the Fermi level in the normal component of the system 
(set Tjano = in eq. I19p . For ((5 < (5p ~ 0.06) however they 
shifts to negative energy opening a pseudogap. Panel D: The 
actual antinodal spectra (Ea„o 7^ in eq. I19p show to become 
asymmetric for 5 > 5c, when the pseudogap starts opening in 
the normal component (see panel C, and bottom right panel 
of Fig. i]). 



of the system (Luttinger theorem^) . The fc„od vector on 
the contrary displays a non-monotonic behavior, showing 
to decrease by reducing doping in correspondence of the 
special doping 5c ^ 0.08, which appears in our study as a 
critical point marking a change in the physical properties 
of the system. In the following section we will show that 
we can relate this behavior of the nodal and antinodal k 
points to a topology change in the Fermi surface. 

In the remaining panels of Fig. [9l we show the spec- 
tral function Im[Gfe(a;)]ii (from eq. [H]). As previously 
explained (see also appendix[X|. within the ED-CDMFT 
method it is possible to display Green's functions on the 
real frequency axis oj by expressing them in a pole ex- 
pansion, displayed by adding in the denominator a small 
imaginary part irj (here we adopted 77 — 0.03i). At the 
nodal point (panel B) the d-wave superconducting gap 
is zero, and a Fermi liquid quasi-particle peak is always 
found at the Fermi level (w = 0) for different doping 5. 
The approach to the Mott insulator {5 0) is marked by 
a progressive reduction of the quasi-particle peaks. This 



behavior is reminiscent of the Mott transition described 
in the standard infinite dimensional Hubbard Model. At 
the antinodal point the superconducting gap is maximal, 
a gap in the spectra is therefore expected. Before looking 
at the full antinodal spectrum, however, it is instructive 
to look at the contribution coming from the normal com- 
ponent, which can provide information on the physical 
properties of the liquid underlying the superconducting 
state. The normal-component spectra can be simply ob- 
tained in our scheme by zeroing the anomalous compo- 
nent of the self-energy 'S,ano{k,uj) in eq. [121 These spec- 
tra are displayed in panel C of Fig. O At doping 6 > 6c, 
a Fermi liquid quasi-particle peak is also found at the 
Fermi level {uj — Q). In this region of the phase diagram, 
the normal properties of the system are therefore Fermi- 
liquid-like (and also knod and kanod are monotonically de- 
creasing with 6, as described above). In correspondence 
of the critical doping 6c however, a pseudogap opens and 
a quasi-particle peak is found at the gap edge at negative 
energy {uj < 0). The pseudogap increases in approaching 
the Mott transition (i5 ^ 0), while, differently from the 
nodal point, the peaks show a roughly constant height. 
This behavior at the anti-nodes of the normal component 
of this superconducting solution can be smoothly con- 
nected to results previously obtained in CDMFT studies 
of the normal state^^'"'^. For 6 < 6c therefore, the nor- 
mal component of the system is not a Fermi liquid in the 
strict sense, at least in the region of momentum space 
close to the antinodal points (but a behavior unusual for 
a Fermi liquid is also detected by the decreasing value of 
knod in panel A) . This behavior appears also in the total 
antinodal spectra (upon restoring the superconducting 
gap Sano in eq. I19p . which we show in panel D. At 
doping 6 > 6c, the quasi-particle peaks present in the nor- 
mal component are parted into two bands by the opening 
of a superconducting gap, resulting in the typical BCS 
symmetric spectra. For 6 < 6c, however, the pseudogap 
already present in the normal component super-impose to 
the superconducting gap, resulting in asymmetric spec- 
tra. This antinodal spectra nicely explain Fig. [6l where 
the local density of states N{uj), directly obtained from 
the cluster-impurity solution, is displayed. The appear- 
ance of the asymmetry for 5 < (5c is therefore interpreted 
by our fc-momentum analysis as the appearance of the 
pseudogap phase, which marks a departure from a Fermi 
liquid based BCS superconductor, once again at the crit- 
ical doping 6c- And these observations must be directly 
linked to experimental spectra, either in the anti-nodes 
of momentum space, obtained for example with angle re- 
solved photo-emissiorii2ii^, and locally in real space, with 
for example scanning tunneling spectroscopy^^. 

In order to characterize the Mott transition, it is useful 
to extract the quasi-particle residuum, which is defined 
as 



dRe^'""'{k,Lu) 
du! 



(33) 



UJ — >0,k—kf 



with kp — knod or kanod in our case. This quantity corre- 
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FIG. 10: (Color online). Quasiparticle residuum at the nodal 
point Znod and antinodal points Zanod as a function of the 
doping 5, evaluated as area of the quasiparticles peaks shown 
in Fig. 1^ (as a crosscheck, Znod is also evaluated as the slope 
of the imaginary part of the self-energy on the Matsubara 
axis, green cross). While Znod monotonically decrease with 
doping, as in the standard picture of the Mott transition, 
the antinodal Zanod shows to stay constant upon opening of 
the pseudogap {5 > 5c ^ 0.08). For comparison's sake, we 
display the nodal Znod obtained by CTQAIC*- on the two 
dimensional Hubbard Model with similar parameters {U = 
12t,t' = 0.0, /3 = 200). These quasiparticle residua should be 
compared with their cluster counter-parts in Fig. |4l 



spends to the area of the quasiparticle peaks (e.g. panel 
A of Fig. [9]), it is unity in the non- interacting case, 
less than unity in an interacting Fermi liquid. In the 
standard description of the Mott transition (i.e. the in- 
finite dimensional Hubbard Model) Z ^ linearly as 
(5^0 (see ref.— ). It is interesting therefore to observe 
how this quantity behave in the nodal and antinodal k- 
points by varying doping S , as we show in Fig. 1101 We 
can give in this way a momentum space interpretation 
of the corresponding cluster Zx that we have discussed 
in Fig. m As we already said, at the nodes well de- 
fined quasiparticle peaks are observed at every doping 5 
at the Fermi level. The residuum Znod is therefore well 
defined according to the expression given above (which 
is strictly valid at w = 0), as in a typical Fermi liquid. 
This is confirmed by the good numerical accord between 
the blue circles (calculated by integrating the area of the 
peaks, which are displayed by introducing the artificial 
broadening parameter 77 = 0.03t) and the green crosses 
(calculated more precise within our method by using eq. 
|33] on the Matsubara axis). For comparison's sake, we 
present also the nodal quasi-particle residuum extracted 
in a CDMFT-study implemented with a different impu- 
rity solver, the CTQMG^, in the two dimensional Hub- 
bard model with Hamiltonian parameters U = 12t and 
t' = Ot. The comparison is only qualitative and it is 
aimed to get insight into the physical trends. We indeed 
observe a monotonically decreasing Znod as a function of 



the doping S, similarly to the standard infinite dimen- 
sional Hubbard Model. We cannot however state within 
our numerical resolution if Znod exactly at the Mott 
point (5 = (as it seems also to suggest the CTQMC re- 
sult), or rather it extrapolates to a finite but very small 
(Znod < 0.02t) value. In the antinodal point more atten- 
tion has to be paid in defining a quasi-particle residuum 
Zanod- As strcsscd above, for S > Sc ^ 0.08 quasipar- 
ticle peaks are present at the Fermi level (panel C of 
Fig. [5]), and Zanod can be well defined by eq. [331 For 
S < 6c ^ 0.08 however a pseudogap opens. Even if a peak 
can be identified at the gap-edge, it is not strictly speak- 
ing a Landau-Fermi liquid quasiparticle, as the imaginary 
part of the self-energy is non-zero (even if small), i.e. the 
quasiparticle has a finite lifetime. Formula [33] cannot 
be directly employed. However we can still calculate the 
area of the peak, and display its behavior as a function of 
the doping S. We find in our result that for 6 < 6c, once 
the pseudogap opens and the peaks move to negative fre- 
quency, the weight Zanod stays constant up to the Mott 
transition point. The behavior of Znod and Zanod here 
presented has to be connected with the effective cluster 
correspondents Zx (Fig. [3]). According to eq. [231 these 
are interpreted in momentum space as the quasiparti- 
cle residua in the corner points of the first quadrant of 
the Brillouin Zone (therefore far from the Fermi surface) , 
while Znod and Zanod, which are instead calculated on the 
Fermi surface, have a real physical meaning. In spite of 
this however, the cluster Zs already embody the physical 
properties (a Z going to zero and another non decreasing 
for (5-^0) characteristic of the Mott transition in this 
two-dimensional system. 

We stress that the description of the Mott transition 
we find in this study of the two dimensional Hubbard 
Model is very different from the standard Mott transi- 
tion picture in infinite dimension^. In our system differ- 
ent regions of momentum space behave very differently in 
approaching the transition point. In this way we can go 
from a Fcrmi-liquid-based superconductor (realized for 
5 > 6c) into the Mott insulator (at 6 — 0) hy passing 
through a phase < 6 < 6p, where the system is at the 
same time insulating in the antinodal region and Fermi 
liquid in the nodal region of momentum space. This lat- 
ter appears to approach the Mott point in the standard 
(infinite dimensional Hubbard Model) way, with a quasi- 
particle residuum Znod (at least within our numeri- 
cal precision) . In the antinodal region instead the quasi- 
particle peak (which underlies the superconducting gap) 
stops reducing a,t 6 = 6c and shifts to negative energies 
opening a pseudogap. This behavior is reminiscent of the 
orbital selective Mott transition, found e.g. in two band 
Hubbard-like models^i^i^S, where the spectral weight is 
not transferred from the low energy (w = 0) to the Hub- 
bard bands (located at a energy scale of order ~ U), 
as in the standard Mott transition, but rather onto a 
smaller energy scale of the order of an exchange coupling 
J ^ l/U'^, inside the Mott gap. In spite our model is a 
one band one, different regions of momentum space ap- 
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FIG. 11; (Color online). Top: Nodal velocity as a function of 
doping measured from the quasiparticle dispersion of different 
materials (figure taken from ref.^^). To be confronted with the 
nodal component of the velocity v^od calculated in our work 
and displayed in the bottom. Bottom: Nodal anomalous ve- 
locity Vano ~ ZT.ano{knod) (taugcut componeut of the nodal 
velocity) and nodal quasiparticle velocity Vnod = ^|Vfc^fc| 
(component of the nodal velocity perpendicular to the Fermi 
surface) as a function of doping 5 (units are aot, where ao 
is the square lattice spacing). The trends are compared 
with CTQMO^ calculations on the two dimensional Hubbard 
Model with similar parameters {U = I2t,t' = 0.0, /3 = 200). 



pear to behave as different bands. By decreasing doping 
5, a first "orbital selective Mott transition" takes place 
at 5c in the antinodal regions, and a full Mott transition 
takes finally place at (5 = 0. 



Useful information on the nodal point can be extract 
by performing a low energy expansion of the Green's 
function (see eq. [21] and [nj) , taking advantage of the 
Fermi liquid properties of the nodal point. In particu- 
lar the nodal velocity has two components, one coming 
from the normal part Vnod ~ kS,k \ parallel to the Fermi 
surface, and the other related to the superconducting gap 
Vnod = ^nod|VfeS;Q„o(fc)|, perpendicular to the Fermi sur- 
face. Vnod can be experimentally extracted, e.g. from 
the ARPES quasiparticle dispersion at the nodeZi, while 
Vano can be determined e.g. as the slope of the super- 
conducting gap at the node^. In Fig. [TT] Vnod and Vano 
are displayed as a function of doping 5. For comparison's 
sake, in the top panel we show experimental nodal ve- 
locity extracted from different materials (the figure has 
been taken from the supplementary material of ref^). 
In the bottom panel we show our result, and also insert 
for a qualitative comparison the CTQMC-CDMFT result 
of refi-1. Vnod shows to be greater than Vano^ in agree- 
ment with experimental observation (see e.g. refJ^). In 
ED-CDMFT it slightly oscillates around a constant value 
from the over-doped to the under-doped [5 < Sc) side of 
the phase diagram, while in the CTQMC case it is slightly 
decreasing with decreasing doping. This behavior is in 
good qualitative agreement with experimental results re- 
ported in the top panel (a quantitative comparison would 
roughly held, for a lattice spacing Oo ~ 4 A and a i ~ 0.5 
eV, Vnod ~ 1 eVA, which is of the order of magnitude 
of experiments). Remarkably, either with ED and with 
CTQMC, Vano displays a dome-like shape, with a maxi- 
mum around optimal doping 6p < S < Sc (notice that in 
the ED case t' = -0.3t Sc 0.08 while in the CTQMC 
case t' — Sc ^ 0.12). This is an important result: as 
shown in formula [30] and discussed in ref4^, the anoma- 
lous velocity Vano = \/2Z„odSano(w = O)smknod can be 
interpreted as a direct measure of the superconducting 
gap in the nodal region, which reveals non-monotonic, in 
agreement with some recent experimental spectroscopy 
results^^i^. The agreement of the trend of Vano between 
the ED and CTQMC solutions shows that this result is 
solid from the theoretical side too. We will come back 
to discuss the nodal and antinodal gap more in detail at 
the end of this section. 

At the nodal point, the combination of the quasipar- 
ticle residuum Znod and nodal velocities can give fur- 
ther information, which can be confronted with experi- 
mental observable quantities and which can further sup- 
port the physical description drawn from our CDMFT 
result. Basing on the Fermi liquid assumption at the 
nodes, in eq. p2p we have for example extracted the low 
energy {ut — > 0) behavior of the local density of states 
N{iL;) ^ ^ — uj. This value is displayed as a function 
of doping § in °Fig. 1121 comparing once again the ED- 
CDMFT results of this work (with t' = -0.3t) with the 
CTQMC-CDMFT results of ref4i (with t' = 0). Once 
again it is the trend we want to compare rather than the 
quantitative values. Starting from the over-doped side, 
the slope of iV(u;) is decreasing monotonically by decreas- 
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FIG. 12: (Color online), 
density of states N{u) ^ 

of doping 5. Trends are compared with CTQMG*" calcula- 
tions on the two dimensional Hubbard Model with similar 
parameters (17 = I2t, t' = 0.0, /3 = 200). 
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ing doping, until showing a slight up-turn close to the 
Mott transition (observable both in the ED and CTQMC 
cases). While the linearity of ^^(0-1) it is well established 
in scanning tunneling experiments^^ , the behavior of the 
slope as a function of doping is at the moment very dif- 
ficult to extract (as it is not possible to obtain absolute 
values for different densities). Our result (in particular 
the up-turn tendency at small doping) has to be therefore 
considered a theoretical prediction. 

A further ratio of experimental relevance can be con- 
nected at first order to the low energy {uj 0) linear 
behavior of the B2g Raman response function XB2g ^ 

- — ^ — UJ (see ref."S^) and the low temperature (T 
0) behavior of the superfluid stiffness Ps{T) — Ps{0) ^ 
— — Lu, which can be extracted from measures of the 
penetration dept h^^i^'^ . Our results (together with the 
CTQMC results of Ref.^^) are shown on the bottom panel 
of Fig. II 31 and compared with the aforementioned Raman 
and penetration depth data presented in the top panel 
(the figure has been taken from refi^). The remarkable 
feature, found in experiments and supported in our calcu- 
lation, is the constant value displayed by this ratio in the 
under-doped region (S < Sc). While in the experimen- 
tal results a sum rule is assumed, in order to being able 
to compare measures from different doping/samples, our 
theoretical results are derived from a bare strongly cor- 
related electron model, where other kind of assumptions 
and approximations (as explained in the previous sec- 
tions) are implied. The convergence of experimental and 
theoretical results therefore strongly supports these find- 
ings, presenting them as distinctive feature of the cuprate 
superconductor nodal dispersion. 

The behavior of the spectra presented in Fig. [§] can 
be directly connected to spectroscopy experiments. In 
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FIG. 13: (Golor online). Top: the linear coefficient of 
the low energy Raman response XB29(<^) ~ Qtj and the 
low temperature superfluid stiffness (from penetration depth) 
ps{T) — ps(0) ~ /3r, extracted from various experiments (fig- 
ure taken from ref.-S^). Quantities are normalized at the op- 
timal doping value. Bottom: The linear coefficient of the 

Raman response and superfluid stiffness (see top panel) as 

^2 

extracted from our calculation — and displayed as a func- 

tion of doping S. Trends are compared with GTQMO^ calcu- 
lations on the two dimensional Hubbard Model with similar 
parameters ((/ = 12t,t' = 0.0, /3 = 200). 



particular, recently a lively debate has risen on the 
momentum resolved structure of the superconducting 
gap^iSi58i59j74^j76^^ We procccd as in ARPES exper- 
iments (see e.gi^), and extract the quasiparticle gap in 
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the nodal and antinodal points of momentum space, tak- 
ing advantage of the periodizing scheme we have intro- 
duced. To this purpose, it is convenient to use the low 
energy (uj — > 0) expansion carried out in eq. 1291 In par- 
ticular we have seen that quasiparticlcs peaks are always 
found in our result, even if not in a strict sense (in the 
pseudogap region 5 < 5c quasiparticle peaks are located 
at the gap-edge and have finite lifetime), and therefore 
we expect the expansion to be reasonably good at small 
frequency (i.e. lo < the superconducting gap). We can in 
this case write the total gap Atot{k,uj) as the quadratic 
sum of two contributions (see eq. 130^ : 



AL (fc, c^) = Al^^ {k,u;) + Al{k,^) 



(34) 



where Asc{k,uj) — Zk^ano{k^uj) is the usual d-wave su- 
perconducting gap (notice that it is directly connected, 
except for constant factors, to the anomalous component 
of the nodal velocity Vano ~ Zk^^^Y^anoij-^ = 0) discussed 
in Fig. [nl), A„o^(fc,w) = |fc(a;) = Zu [6 +ReI]fe(w)] is 
a normal contribution to the gap which can arise only if 
for some k and w — > 0, the normal component of the self- 
energy grows enough so that the band equation 
S,k — Q cannot be satisfied, i.e. there is not Fermi sur- 
face. Now, this does not take place at the nodes, where 
quasiparticlcs are found for all dopings (panel B of Fig. 
[9]), ^fc = 0, the system presents standard Fermi liquid 
properties, and the total spectral gap coincides with the 
superconducting gap Atot = Age- (This already suggests 
that in looking for the "real" superconducting gap of the 
system one should look at the nodal region gap. Re- 
cent theoretical^^ and experimental Raman spectroscopy 
studies^ point at this direction.) In the antinodal re- 
gion, instead, we have observed that a pseudogap opens 
in the normal component for 6 < 0.08 (panel C of Fig. [9]), 
and this fact is associated in our calculation with the ap- 
pearance of lines in fc-space close to the antinodal regioE 
where the self-energy is diverging (see refi2S.i4^). It is not 
possible therefore to satisfy the equation = 0, and a 
Fermi surface does not exists anymore. The normal con- 
tribution Anor kicks in, and determines the properties oi 
the total antinodal gap Atot , originating the asymmetric 
spectra we already described in panel D of Fig. [SI 

We can extract, like in ARPES experiments (e.g. 
ref.-^), the antinodal gap Atot from the spectra of panel 
D of Fig. [5] by measuring the distance of the quasipar- 
ticle peaks at the gap-edge from the Fermi level uj — 0. 
In the same way, from panel C of Fig. [9l we can mea- 
sure the normal contribution Anor, and display them as 
a function of doping S. We can compare directly with 
experimental results, which are shown in the top panel 
of Fig. [M] (for convenience's sake we extract the picture 
from the Raman results of ref.^^, but ARPES points are 
also displayed), while in the bottom panel we show the 
results of our calculation. The antinodal gap Atot is a 
monotonic decreasing function of doping (curve labeled 
"Anti Nodal Big" in the top panel of Fig. [H)) . as it 
has been known from experiments since a long time (see 
e.g. ref^i^) and predicted in the most popular theories 
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FIG. 14; (Color online). Top: Nodal (Big) and antinodal 
(B29) quasiparticle gap extracted by different spectroscopy 
experiments (in particular, the figure is taken from the Ra- 
man spectroscopy results of refi^) as a function of doping. 
Bottom: The antinodal gap Atot and the nodal Asc ex- 
tracted from the spectra of Fig. [9l as a function of doping S. 
Anor is the pseudogap in the normal component extracted in 
panel C of Fig. H 



of high temperature superconductivity ( e.g. resonat- 
ing valence bond theory"*, for a recent general review 
see e.g4^). In our calculation we show that at small 
doping an important contribution comes indeed from the 
normal component Anor, which appears at Sc and it is 
also monotonic. From the experimental side, the nov- 
elty comes from precise measures of the nodal gap, re- 



cently obtained by Raman spectroscopy^^ and ARPES^^, 
which show surprisingly that in this region of momentum 
space the spectral gap is non-monotonic with doping S, 
tracking instead the behavior of the critical temperature 
Tc (curve labeled "Nodal 623" in the top panel of Fig. 
[T4| . In our calculation the nodal gap corresponds to the 
anomalous component of the nodal velocity Vano, which 
we have discussed in Fig. 1111 and which shows indeed a 
behavior strongly similar to these experimental results. 
We clarify now how it is possible that the nodal com- 
ponent of the gap tracks Tc as a function of doping S, 
while at the same time the antinodal component is mono- 
tonic, by disentangling he superconducting contribution 
Asc = \/ Ajot ^ ^nor; which wc also display in the bot- 
tom panel of Fig. [TH To check the validity of our formula 
1341 and making connection with Vano ~ Zk^^^'S'ano, which 
has been evaluated at the nodal point, we also evaluate 
and display A^c = Zanod£'ano{kanod,t^ = Anor), finding 
numerical agreement. This shows indeed that the total 
antinodal gap Atot has indeed two distinct contributions, 
displaying opposite trends with dopings. In the under- 
doped side A„or dominates in the antinodes, and creates 
a monotonically increasing total gap. In the nodes, in- 
stead, A„or is zero and only the superconducting gap is 
detectable. In our results we connect therefore the exper- 
imentally observed two-gap phenomenon with the open- 
ing of the pseudogap at the anti-nodes at a finite critical 
doping (5c. And, according to our point of view, this is 
an effect arising in a strongly correlated electron system 
that approaches the doping-driven Mott transition in a 
two dimensional lattice. The fact that A^c tracks T^S^ 
is remarkably similar to the standard BCS superconduc- 
tivity. These results put strong constrains on the the- 
ories of cuprate-based superconductivity, which have to 
consider the presence of these two distinct components 
in the spectral gap and its interplay with the dome-like 
shape of the order parameter, the rising of the pseudogap 
phase and the approaching to the Mott insulator. Our 
cluster results, interpreted via a periodization procedure, 
well fit the experimentally observations and give a sim- 
ple interpretation in terms of a combination of all these 
effects. 



V. A MIXED-PERIODIZATION SCHEME 

The discussion presented in the previous subsections 
is valid close to the nodal and antinodal points of mo- 
mentum space, where we have shown that periodizing 
the self-energy or the cumulant its a reasonable approxi- 
mation. Obtaining information in intermediate region of 
momentum space, between the nodal and antinodal ones, 
is beyond the limits of a pure 2 x 2-plaquette study. In 
this region in fact the Fermi liquid physical properties of 
the nodal points have to interlace in some non-trivial way 
with the insulating-like properties of the antinodal point. 
A detailed description of this phenomenon can be taken 
into account by studying bigger cluster (i.e. obtaining a 
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FIG. 15: Left: the shaded region is the patch $fc used in the 
phenomenological Boltzmann approach of ref.~ in the first 
quadrant of the Brillouin Zone. Right: the patch used in 
this work (see formula |37|I. 



better truncation of the Fourier expansion in eq. [20)) . at 
present not possible with ED-CDMFT. Nevertheless we 
want to keep in this section a low profile, and introduce 
a first order description of the spectral properties in all 
momentum space, which could be compared to experi- 
mental results (like e.g. ARPES). We want therefore to 
introduce a periodizing scheme able to describe 

1. a Fermi liquid quasiparticle in correspondence of 
the nodal point 

2. the opening of a pseudogap in the antinodal points 
in the under-doped regime 

3. the formation of the Mott gap at high energies in 
approaching the Mott insulator (see top panel of 
Fig. [8] region marked A). 

We can satisfy the first condition by using the self-energy 
S-periodization. The second and third points are instead 
obtained by using the cumulant A^-periodization, which 
is able to describe the formation of the antinodal. 

In order to combine all these requests, we take 
advantage from ideas introduced in phenomenologi- 
cal Fermi-liquid Boltzmann approaches to the normal- 
state transport properties of cuprate superconducting 
materials^i^i^i^. We base in particular on the work 
of ref4^, where the division of momentum space in 
nodal Fermi-liquid-like regions, cold spots, and antinodal 
insulator-like regions, hot spots, was achieved by project- 
ing the quasiparticle scattering operator Ckk' on a basis 
of patches {(j)a\ in momentum space with different tem- 
perature scattering-dependencies: 



kk' 



J2 Mk)Cc.p{T)<j}p{k) (35) 

a/3 



In this case, it is possible to solve exactly the Boltzmann 
equation for the simple case of two patches (one marking 
the nodal and the other the antinodal region). The shape 
and the scattering properties of the cold patch (modeled 
by a small set of parameters) was fixed by obtaining the 
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best fit on few transport quantities (resistivity, Hall co- 
efficient), and a systematic correspondence with others 
transport quantities (magnetoresistance and termoelec- 
tric power) was then obtained. For comparison, the nodal 
patch used in ref."'^ is shown in the left hand side of 
Fig. [15] (the antinodal patch is simply defined as 1 — 

In the same spirit, we introduce here a mixed 
periodization-scheme, by projecting the lattice self- 
energy Sfc(ti;) on a nodal <i>(fc,w) and an antinodal 
1 — $(fc,Li;) patch in momentum space: 

Efe(w)= $(fc,w)Sfe(w)[S] + [l-$(fc,w)]Efe(u;)[>l] (36) 

where I]fe(cL;)[S] is obtained by periodization of the clus- 
ter S and Efe((jj)[A^] by periodization of the cluster cu- 
mulant M (eq. [^5)1 . The patch ^{k,uj) separates the 
E-periodized from the A^-periodized regions (right hand 
side of Fig. [T5)) . For convenience's sake, we choose a form 
with k and u separable: 



^{k,uj) = *(fc)«'(w) 
*(a) = ^ [1 - tanh(/3„r„)] 



(37) 



where is the distance from the a center in fc or w space: 

rk = \J [kx - fcrj^ + (ky - ky^Y ~ ro{9) 

= \ll>\ - ll>o (38) 



Here (fc.,,fc,J ~ (f,f), r„(0) 



To e 



and e = 



a,Tcta,n{ky / kx) can be chosen to properly shape the patch. 
We fix through all the rest of the paper ujo ~ 0.5i, f3k ^ 
100 and ~ 15, which smooth the boundaries of the 
patch. It remains to fix the parameters kc, ro and CTq, 
which have to be chosen to mimic the phenomenological 
patch of Fig. [151 There is of course a good degree of 
arbitrariness in its form. We want however to have a first 
order qualitative description of the physical properties in 
momentum space, with the aim to address consideration 
that are only slightly dependent on the exact form of the 
patch. 



A. Evolution of the Fermi Surface 

In order to fix the dimension of the patch 5'(fc) in k- 
space it is useful to follow the evolution of the Fermi 
surface in approaching the Mott insulator with the two 
periodization schemes E and M, as shown in the panels 
of Fig. [ini At high doping (in our case 6 = 0.13) all 
the system is well described by a Fermi liquid, the self- 
energy is mostly local, the E— (green dash line on the 
left top panel of Fig. [TB)) and M— (blue continues line 
on the left top panel of Fig. [TH]) periodizations give in 
practice the same result. We can choose to describe the 
system with the E-periodization and ^'(fc) covering all 
the fc-space (for example see the left panel of Fig. [TTlfor 
6 = 0.13, where the patch is marked by the gray region 
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FIG. 16: (Color online) . We show the evolution of the normal- 
component Fermi surface (blue continuous line in the Ai- 
periodization scheme, green dashed line in the E-periodization 
scheme), defined has + Re Efc = in eq. 1191 Panels span 
from the over-doped region {S — 0.13 top- left panel) to the 
Mott insulator (S — 0.0 bottom-right panel). By periodiz- 
ing the cumulant M, the Fermi surface disappears in the 
insulating state and it is replaced by a line of zeroes of Gt 
(red lines). The way this takes place is via a topological phase 
transition around Sp ~ 0.06 < 5 < <5c ~ 0.08, from a large 
hole-like Fermi surface to a pocket-like Fermi surface. The 
low doping region is marked by a phase where Fermi surface 
and lines of zeroes co-exists. This gives origin to the pseudo- 
gap and the fragmentation of the Fermi arcs in measures of 
spectra of the normal state system (see e.g. Fig. I18|l . By pe- 
riodizing the self-energy E instead, by construction, the 
lines of zeroes cannot appear. The pseudogap is realized only 
trough a modulation of the spectral intensity along the Fermi 
surface line, which extends always on all the Brillouin Zone 
(see also Fig. [T5)) . 



covering all the quadrant). By reducing doping, however, 
the off-diagonal components of the cluster self-energy S 
are not negligible anymore, and the two periodizations 
produce different results. In particular, as we stressed in 
the previous sections, doping Sc ^ 0.08 and 5p ^ 0.06 are 
special points. In the Al— periodization the Fermi sur- 
face shows a striking topological phase transition, pro- 
duced by the appearance of lines of zeroes of the Green's 
function Gk {(^ — > 0) at the Fermi level (marked by a con- 
tinuous red line in Fig. [TB)) . The effect of the appearance 
of the lines of zeroes of the Green's function is at the ori- 
gin of the opening of a pseudogap in the spectral function 
close to the (0, tt) — (tt, it) and (tt, 0) — (tt, tt) sides of the 
first quadrant in the Brillouin zone (as shown for exam- 
ple in panel C of Fig. [9]) . In the E-periodization instead 
a more continuous evolution of the original high-doping 
Fermi surface (green dashed line) takes place at all dop- 
ings up to the Mott insulating state, where it disappears. 
We notice however that the doping 6c still marks a change 
in the curvature of the Fermi surface (even if with the E 
periodization the effect is more difhcult to be noticed at 
naked eye) , as evidenced in the downturn of the knod vec- 
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tor as a function of doping in panel A of Fig. [9l This 
effect was first noticed in refi^, where we show that the 
resuh of periodizing E produces a Fermi surface which 
enhances its hole-hke curvature while reducing doping 
(and the spectral weight reduces too with respect to the 
nodal point). This goes in the direction of forming a hole 
pocket, which however never arrives to be created within 
the E— periodization. Within the M scheme instead, the 
Fermi pocket forms at low doping and its progressive re- 
duction in approaching the Mott transition describes the 
way the Fermi surface disappears. In the Mott insulat- 
ing state Fermi lines have of course totally disappeared, 
but a line of zeroes of Gk (i-e. a line of divergent self- 
energy S^°''(a; = 0)) remains (in red) close to the corner 
k = (tt, tt) of the first quadrant of the Brillouin Zone. 

It is very difficult to state how close to the real solution 
one or the other of the two descriptions are. Particular 
intriguing is the M. periodization result. Not only it well 
portrays the pseudogap in anti-nodes, it also produces, 
together with the Fermi surface, lines of zeroes of Gk ■ In 
this way it describes a continuity from the Fermi liquid 
at high doping (where only the Fermi surface is present) 
to the Mott insulating state (where only lines of zeroes 
are present). To this respect, it is clear that the E pe- 
riodization scheme fails in describing the Mott state, as 
it is unable by construction to build up lines of zeroes. 
Within our analysis however, we cannot claim that in 
the real physical system a Fermi pocket, together with 
lines of zeros, is actually present. Recent experiments on 
cuprate systems, where it has been possible to induce a 
low temperature normal state by the application of an 
external magnetic field, have actually observed de Haas- 
van Alphen oscillations compatible with a Fermi pocket 
picturo^^. Hall resistivity measures, extracted at low 
temperature by suppressing superconductivity with the 
application of an external magnetic fiel d^°'^^i^^ , are also 
compatible with the scenario of a topological phase tran- 
sition of the Fermi surface. Our results are also in strong 
resemblance with the theoretical study of refi^, where 
similar conclusions on the evolution of the Fermi surface 
and the appearance of lines of zeroes have been drown 
starting from an ad hoc model for the doping-driven 
Mott insulator transition in two dimensions. Other non- 
perturbative microscopic approaches have drawn conclu- 
sions in similar directions (see e.g. refi^i^). There 
are however some caveats in concerning the hole pocket 
which have to be considered. 

In first place, according to the generalized Luttinger 
theorem^i^iS^iSl the volume enclosed between the Fermi 
surface and the line of zeroes (if present) should be equal 
to the particle density n — 1 — 5. This theorem is quite re- 
spected (but not so much at low doping) by the Fermi sur- 
face derived with the E-periodization (green dashed line 
of Fig. [in]) . It is not clearly obeyed instead at low dop- 
ing by A^-periodization, as evident e.g. by looking at the 
area enclosed between the Fermi pocket (blue line) and 
the lines of zeroes (red line) in the panel 5 = 0.05 of Fig. 
[TBI (close to half- filling the volume which gives the correct 
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FIG. 17: (Color online). The dimension of the patch 'ilk in 
momentum space (gray shaded region) is determined in such 
a way that the A^-Fermi surface (black line) is cut away and 
replaced by the E-Fermi surface (dashed green line). The red 
line are where the self-energy is diverging and Gk = 0. 



density should be half of the quadrant). Far from stating 
that the system is violating Luttinger theorem (which in 
strongly correlated system is actually a possibility's,;^), 
this effect is likely an artifact coming from the truncated 
Fourier expansion in eq. I20| which can be improved only 
by increasing the cluster size. A possible scenario is that 
in real systems the line of zeroes is closer to the zone 
diagonal which goes form (0, tt) — > (tt, 0), as actually we 
find at half-filling for a chemical potential value in the 
middle of the Mott gap /j, ~ C//2 (the panel displayed in 
Fig. [in] has a chemical potential close to the Mott gap 
edge /i ~ J7/2 — Am/2, where Am is the total Mott gap). 
This is in fact what was propose in ref.— as starting hy- 
pothesis. In this way, the side of the hole-pocket facing 
the (tt, tt) corner point is "cancelled" by the proximity of 
the lines of zeroes (see for example the spectral functions 
of Fig. I17p . and the resulting picture is a Fermi arc in 
the nodal region, which is replaced by a lines of zeroes in 
the antinodal regions. 

In second place, the discussion we have previously car- 
ried out (see Fig. [5|in section [Til] and the description of 
nodal properties in section [TV]) shows that the low energy 
nodal point is better portrayed by periodizing S, i.e. at 
the nodal point we rather have a Fermi arc more than a 
pocket. This is important if the nodal point properties 
(presented in the previous sections JTlT] and ]TV]) have to 
be well portrayed. The real solution result at the nodes 
is likely to lay in between the M and E periodization 
schemes. 



B. The choice of the patch shape 

In line with the discussion above, we choose therefore 
to assign to the nodal region the path ^E'fe, which uses 
the E periodization. The rest of the A:-space (covered 
by 1 — '^k) is described with the A1 periodization. By 
following this criterion we establish the evolution of '^{k) 
with doping, as shown in Fig. [TT] The size and shape of 
the patch '^k is determined so that the M Fermi surface 
is cut away, and the only piece of Fermi surface inside 
^'fe is the one produced by periodizing E (green dashed 
line). By reducing doping 5 the size of the patch '^k 
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FIG. 18: (Color online). The spectral function — ^ImGfe(LJ —>■ 
0) in the first quadrant of the Brillouin Zone obtained with the 
mixed periodization (top row) is confronted with one obtained 
from the self-energy (E) and the cumulant (Ai) periodizations 
(second and third row from the top). The broadening rj — 
0.052, the color scale x = 0.5 (see appendix El for details). 

progressively reduces, until possibly disappearing in the 
Mott insulating state {6 = 0) where only the Al-scheme 
reproduces well the gapped spectra. We have fixed the 
center of ^ffe by choosing {K^,kyJ = (^a;„odi ^y^od); the 
momentum coordinate of the nodal point described in 
panel A of Fig. [HI A complete description of the patch 
parameters Vo and Uo (see formula 1571) as a function of 
doping is given in the following table: 
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We notice that with the choice of this patch most of 
the lines of zeroes disappear from the quadrant the Fermi 
level. As mentioned above, we are not able to make 
definitive statement about the actual position of these 
lines of zeroes. If it were closer to the quadrant diagonal 
(as conjecture in work—), they would re-appear in the 
1 — ^'fc region of momentum space, which is described 
by the cumulant. As stressed in the previous subsection, 
this would be important if one want to respect the Lut- 
tinger theorem on the particle counting. In this mixed 
periodization used for the 2x2 plaquette result however 
this does not take place. We remark however that in the 
1 — '^k region scattering rate (i.e. self-energy) is very 
high (similarly to the results in refi^), and this fact is 
ultimately the reason for the appearance of a pseudogap. 





FIG. 19: (Color online). Comparison between the E- 
periodization scheme, the cumulant A^-periodization scheme 
and the mixed-periodization scheme for a low doped system 
5 = 0.05. On the top row we show the full energy range 
— 14 < < 14 covering the lower Hubbard and upper Hub- 
bard bands. The gross features are very similar. Notice the 
E-scheme introduces fake density of states in the Mott Hub- 
bard gap. In the bottom row we show a close up at low 
energy — 1 < u < 2, where the methods are most different. 
The mixed scheme is constructed to retain a E-periodization 
character close to the nodal region k ~ (f , f ) and a M- 
periodization character in the antinodal region k ~ (0, tt). 
Spectral functions are displayed by introducing a uj depen- 
dent broadening irjiuj) and the maximal color scale value are 
X = 0.20(0.25] for the top [bottom] panel (see appendix IX)!. 



whether or not the lines zeroes are effectively present. 
The division of the fc-space via the patch ^I^fc is surely 
artificial, and it has not pretext of describing in detail 
the real properties of the system. It has however capa- 
ble of capturing in a unique scheme either the virtues 
of the S-periodization (above all in the nodes) and of 
the Al-periodization (above all in the anti- nodes), which 
portray with some good degree of confidence the physical 
properties in different regions of momentum space (as we 
discussed in the previous sections). 

A confront between the E, M and mixed periodiza- 
tions is presented in Fig. ITBland Fig. [111 In Fig- HHwe 
show the spectral density A{k,Lo) = ilmG^°''(w — > 0) 
in the first quadrant of the Brillouin Zone for decreasing 
doping 5 (from left to right). The mixed-periodization 
scheme (top row) is confronted with the S and pe- 
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riodizations (bottom rows), showing how the patch 
is interlacing them. The well known phenomenon of the 
Fermi arc breakup is reasonably well described by all 
methods, showing this is a solid result of CDMFT. More- 
over, the similarity with spectra calculated in previous 
CDMFT work on the normal state^"^ shows the smooth 
continuity between the normal component spectra of this 
superconducting state solution with the spectra of a nor- 
mal state result. In Fig. [19] we show the band spec- 
trum A{k, uj) plotted as a function of the energy co along 
the path (0,0) (7r,7r) -> (0, tt) -> (0,0) in the first 
quadrant of the Brillouin Zone. In the top panel we con- 
front the full energy range — 14 < w < 14 covering upper 
Hubbard band and lower Hubbard band. At this energy 
resolution the three schemes are qualitatively very sim- 
ilar. We just stress that the S-periodization artificially 
introduces spectral weight in the Mott gap (as discussed 
in Fig. [8] and evident in the figure around u ^ 3 — It). 
This justify the choice of cutting the patch '^{lo) at low 
energies uj < 0.5t, using the A^-periodization in the re- 
maining of the energy range. In the bottom panel of 
Fig. lfTHl) we show a close up at low energy (— 1 < < 2). 
In this energy range the methods most differ, however the 
qualitative results are still very similar. In particular the 
mixed scheme has been designed to well describe a Fermi 
liquid linear dispersion at the node, and the right descrip- 
tion of the Mott gap, especially in the antinodal region. 
In the analysis of the following sections we will therefore 
apply the mixed scheme introduced here through formula 
1361 to periodize the normal component of the self-energy 
(while the anomalous component Sa„o(w) is always ob- 
tained through formula I24i implying by construction a 
d-wave shape of the superconducting gap). 



C. Local density of state with the mixed scheme 

We use now the mixed periodization to re-calculate 
the local density of states N{uj) = —-^ J2k IniG'fc(ijj) at 
low energy, beyond the cluster energy resolution (which 
is due to the finite dimension of the truncated Ander- 
son impurity model used to implement CDMFT, see Fig. 
E]) . As widely explained above, this relies on the implicit 
assumption that the superconducting gap has a d-wave 
form (see eq. [24]) and at the node we have a well de- 
fined Fermi liquid arc. The result is shown in Fig. 1201 
A "V-shaped" N{uj) is observed either in the over-doped 
and under-doped regions. Coming from the over-doped 
side towards the under-doped side, the slope is always 
decreasing, until reaching a saturating value at the small 
doping (until eventually showing a small up-turn for the 
smaller doping, see Fig. [T^ . The slopes well fit the ana- 
lytical value extracted at low energy in formula 1321 On a 
wide range of energy (— 1 < a; < 1), the V-shape is quite 
symmetric in the over-doped region 6 > 6c ^ 0.08 (as al- 
ready remarked in the CTQMC-CDMFT study of ref4l). 
Spectra become strongly asymmetric in the under-doped 
region(5 < Sp 0.06), when the pseudogap in the nor- 




FIG. 20: (Color online). Local density of states N{uj) as a 
function of the doping obtained with the mixed-periodization 
scheme. The V-shape at w — > is expected from a d-wave 
superconductor with nodal Fermi liquid quasiparticles. For 
doping S < Sp 0.6 spectra develop a marked asymmetry. 
This result interprets the cluster result of Fig. [6] and it is in 
good agreement with scanning tunneling experiments (see e.g. 
ref.— ). The display on the real axis of the Green's function 
has been obtained by introducing a tj-dependent broadening 
ir;(w) (see appendix \K\ for details). 



mal part of the system opens and super-impose to the 
superconducting gap (as commented in panel D of Fig. 
[H). The qualitative behavior of these curves should be 
compared with the raw cluster result of Fig. [B] The 
periodizing mixed scheme we introduced should be con- 
sidered as the best fit we could achieve to our cluster 
DMFT results, basing on few reasonable solid assump- 
tions on the physical properties of the system (like the 
d-wave superconducting gap, the pseudogap formation, 
the Fermi liquid properties at the node). This allows us 
to recover a momentum dependent Green's function and 
access physical quantities comparable to experimental re- 
sults on cuprates. The qualitative behaviour of N{ll!) we 
determined well portrays in fact results of scanning tun- 
neling experimental^, and it supports a comparison in 
momentum space. 



VI. QUASIPARTICLE SPECTRA FROM THE 
SUPERCONDUCTOR TO THE INSULATOR. 

In this section we use the mixed-periodization scheme 
and derive a detailed description of the connection be- 
tween the spectra of the Mott insulator and the super- 
conductor by varying doping. 
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FIG. 21: (Color online). Band dispersion A{k,uj) along the 
path (0,0) — > (7r,7r) (0, ti") (0,0) in the first quadrant 
of the Brillouin Zone. The Mott-insulating state (5 = 0) is 
confronted with a small doping {5 — 0.04) state. Top-left 
corner: the full range of energies ( — 15 < cj < 15) covered 
by the one-particle spectrum is displayed. The gross band- 
structure, displaying a lower Hubbard band (LHB) (—6 < 
ui < 0) and a upper band (UHB) (10 < a; < 15) separated by 
a Mott gap (0 < a; < 10) remains un-changed in going from 
the insulator into the metal. Top-right corner: detail of the 
LHB, which rigidly shifts in going from the insulating to the 
metallic state. Bottom-left corner: detail of the low energy 
(— 1 < tj < 2) feature. In going into the metallic state, the 
shift of the bands is not rigid and spectral weight appears at 
positive energy as soon as doping is added into the system. 
Bottom-right corner: detail of the UHB. Upon doping, a 
rigid shift in the metallic state is accompanied by a transfer 
of spectral weight to low energies. The color scale maximum 
value is always x — 0.20, except for the lower left panel, where 
it is set X = 0.25. We refer to appendix 1X1 for the choice of 
the broadening irj{Lu). 



A. Doping a Mott insulator 

To this purpose, we analyze in Fig. [21] the spectral 
functions, A{k,uj) = —-lmG{k,uj) in the uj — k space, 
along the fc-path (0,0) ''-^ (7r,7r) (0, tt) (0,0) in 
the first quadrant of the Brillouin zone. We compare 
the insulating state {5 = 0) with a slightly doped state 
[5 = 0.04). In the left top panel we display a wide 
energy-range — 14 < < 14, which covers the lower 
(LHB) and upper (UHB) Hubbard bands. We notice 



that, upon adding a small doping, the gross structure of 
the lower Hubbard band (—6 < w < 0) and the upper 
band (10 < a; < 15), which are separated by a Mott 
gap (0 < < 10), remains substantially un-changed. 
A detail of the lower Hubbard band, a closeup of the 
Fermi level (— 1 < w < 2) and a detail of the upper 
Hubbard band are presented in the right top, left bot- 
tom, right bottom panels respectively. The LHB (top 
right) rigidly shifts in going from the insulating to the 
metallic state, loosing spectral intensity, which goes to 
build up quasiparticles at the Fermi level ld = Q (bottom 
left). The novel outcome from the CDMFT calculation, 
as compared with the sing le-site DMFTi^, consists in the 
anisotropic fashion quasiparticles first occupy the Fermi 
level in momentum space. This is better seen in a closeup 
of the band dispersion around the Fermi level (bottom 
left). The region where firstly quasiparticles appear is 
close to A; = (^, -1). In fact, already in the insulating 
state {6 — Q), we observe that around k — (^, ^) a heavy- 
particle hook-shaped band is closest to the u ~ Q level, 
while in the proximity of A: = (0,7r) there is a "pseudo- 
gap" and the band disperses at negative energies. This 
kind of band-structure survives upon doping, as we can 
observe in the 5 = 0.04 panel. Contrary to the Hub- 
bard bands, however, the low-energy band does not shift 
rigidly as doping is added to the insulator, rather it stays 
pinned at the Fermi level, and the shifts in frequency is 
only a small fraction of the changing in chemical poten- 
tial. We also observe the appearance of spectral weight 
at positive energy, coming from both the LHB and UHB, 
which starts building up a full Fermi-liquid-like band, as 
we show more in detail in the following. In the bottom 
right corner of Fig. [21] we finally show the UHB. In this 
case a rigid shift in the metallic state is accompanied by 
a strong reduction of spectral weight. In fact, accord- 
ing to our result, the UHB narrows with respect to the 
insulating state. 



B. Approaching the Mott insulator from the 
over-doped side 

We start now from the viewpoint of the highly doped 
system and observe how the approach to the Mott insu- 
lator affects the electronic structure of the superconduct- 
ing state. The high-doping system offers the advantage 
of having more standard Fermi-liquid-like properties (the 
patch ^'(fc) introduced in the previous section covers all 
the momentum space) . In Fig. [22] we show the elec- 
tronic band A{k,u}) in the k — uj space, once again in 
the path (0,0) (tt, tt) -> (0,7r) -> (0,0) of the first 
quadrant of the Brillouin zone, on the full energy-range 
10 < w < 15 covering the LHB and UHB. Here we want 
to display the evolution from high doping (left) to small 
doping (right). At high a doping 6 = 0.13 a Mott Hub- 
bard gap which separates a LHB structure from the UHB 
is already visible. A narrow but Fermi Liquid-like band 
is however present, and it crosses the Fermi level either 
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FIG. 22: (Color online). Evolution of the one-particle spec- 
trum from the over-doped (left side, 5 = 0.13) to the under- 
doped (right side, 5 = 0.04) state. The typical Fermi-liquid 
band, well visible in a narrow energy range (— l<a;<l) 
around the Fermi level at high doping S = 0.13, is progres- 
sively destroyed by reducing doping, while the upper and 
lower Hubbard bands build up acquiring spectral weight . 
These features are shown in detail in Figures [23l and [24l The 
color scale maximum value is x = 0.20 and we refer to ap- 
pendix |X] for the choice of the broadening iTj{uj). 



in the region of momentum-space around = ^) and 
k = (0,7r) (see also Fig. [24]). Reducing the doping 6, 
we see that this narrow band loses intensity more and 
more (follow the horizontal line at fc = (tt, tt) in o; ^ 0) 
to the advantage of the Hubbard bands, which instead 
gain spectral weight in approaching the Mott insulator. 
To have a glance on how this is taking place we look in 
the following figures at the different energy-regions of the 
band in further detail. 

Fig. m shows the UHB and the LHB in detail. The 
statements made above are confirmed: both Hubbard 
bands gain spectral weight in decreasing doping, mainly 
in the region of momentum space close to the anti-nodes 
k ^ (0,7r) (follow once again the horizontal line), while 
rigidly shifting with respect to the change of chemical 
potential A/i (one can actually show that the shift of 
the bands in energy is equal to A/i). The behavior of 
the Hubbard bands is therefore in agreement with a pic- 
ture describing the approach to the Mott insulator as a 
rigidly moving bands, which transfer part of their weight 
to low energy. As mentioned in the previous subsection, 
the novelty of our CDMFT result in finite dimension, 
with respect to the standard vision of the Mott transi- 
tion given by single-site DMFT in infinite dimension, is 
that this transfer of spectral weight takes place in a very 
anisotropic fashion, with the antinodal regions k ^ (0, tt) 
getting insulating before the nodal ones. 

A closer look to the behavior at low energy. Fig. [Ml 
reveals much richer phenomena taking place. In the top 
panel we show the spectra resulting from the one-particle 
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FIG. 23: (Color online). Top panel: Evolution of the LHB 
(top) from the over-doped (left side, S = 0.13) to the under- 
doped (right side, 5 = 0.04) state. The LHB rigidly shift 
according to the change in chemical potential by reducing 
doping, acquiring structure (follow the horizontal line at fc ~ 
(0, tt). Bottom panel: Evolution of the UHB from the over- 
doped (left side, 5 = 0.13) to the under-doped (right side, S = 
0.04) state. The UHB rigidly shift according to the change 
in chemical potential by reducing doping and, as evident in 
the color intensity, it gains spectral weight. The color scale 
maximum value is a; = 0.20 and we refer to appendix [XI for 
the choice of the broadening iri{uj). 



Green's function Gii{k,uj) (eq. [TO]) in the superconduct- 
ing state. First of all, contrary to the behaviour of the 
Hubbard bands, the low energy part of the band does 
not shift proportionally to A/i, in agreement with the 
observation already made in Fig. 1211 The Fermi-liquid- 
like band at high doping (left side, 6 — 0.13) is progres- 
sively destroyed by decreasing doping (with progressive 
reduction of spectral weight in the arc at < w < 2). 
In the region close to k — (0, tt) instead, the d-wave 
superconducting-state opens a gap by removing spectral 
weight from positive and negative energies (a Bogoliubov 
band is formed at > 0). As widely discussed in the 
previous sections however, the presence of lines of zeroes 
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FIG. 24: (Color online). Top panel: we show the evolution of 
the low-energy one-particle spectral band from the over-doped 
(left side, S = 0.13) to the under-doped (right side, S — 0.04) 
state. A Fermi-liquid-like band at high doping (left side, S — 
0.13) is progressively destroyed by decreasing doping. Around 
the region k = (f , -f) the dispersion presents a Fermi liquid 
quasiparticle crossing the Fermi Level. Around the region 
k — (0,7r), the d-wave superconducting- state opens a gap by 
removing spectral weight from the Fermi level to positive and 
negative energies. Bottom panel: we show the contribution 
to the low energy band coming from the normal component of 
the system (we set the anomalous self-energy Eano = 0). In 
particular we observe that in the under-doped side {5 < 0.08) 
a gap (the pseudogap) opens in the k — (0, tt) region even if 
the superconductive term is absent. For convenience's sake, 
the color scale maximum value is set x — 0.30 for S = 0.13 
and X — 0.25 for all the other doping values. We refer to 
appendix \K\ for the choice of the broadening ir;(tj). 



in the Green's function kicks in the under-doped region 
((5 < 0.08), opening a pseudogap in the normal compo- 
nent of Gii(fc,tt') around k — (0,7r). This is shown in 
the bottom row of Fig. [Ml where the same panels of the 
top row are reproduced by imposing the anomalous self- 
energy identically zero Ttano = in eq. \TU\ Switching off 
superconductivity has in general little effect on the quasi- 
particle bands, except indeed close to the antinodal point 
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FIG. 25: (Color online). The the low energy kink-feature in 
the spectral intensity at the nodal point for various values 
of doping. For convenience's sake, the color scale maximum 
value is set x = 0.35 for 5 — 0.13 and x = 0.30 for all the 
other doping values. We refer to appendix |X] for the choice of 
the broadening irj{Lo). 

k ^ (0, tt) and close to the Fermi level w ~ (the nodal 
points fc ~ (f , f )) are practically unaffected). The differ- 
ences are most evident in the over-doped region (see e.g. 
6 ~ 0.13 — 0.08 close to /c ~ (0,7r)), where the supercon- 
ducting gap has disappeared and the quasiparticle band 
is reconstructed. The weight in the Bogoliubov branch 
at a; > also disappears. For S > Sc ^ 0.8 differences 
are much less evident (above all S — 0.04). A gap is also 
present in the normal component solution, the structure 
of the quasiparticle dispersion around k ^ (0, tt) is very 
reminiscent of the superconducting solution in the cor- 
responding top row. It seems that in the under-doped 
region the superconducting gap appears to complete the 
structure already present in the normal component so- 
lution. This indicates that the pseudo-gap appearing 
in the normal component of the system (and present at 
temperatures above Tc) and the d-wave superconduct- 
ing gap are possible answers of the system to the same 
instability, and they co-exist in the under-doped region 
(while in the over-doped only the superconducting gap 
is present). This instability is, in our view, connected 
to the approach to the Mott transition, and it reflects 
the anisotropic way chosen by the Hubbard system in 
two dimension to approach the insulating state: the re- 
gions around the anti-nodes k ~ (0, tt) become insulting 
(at 6 = Sc ^ 0.08) before the regions close to the nodes 
k (f I f ) (which become finally insulating (5 = 0). 

VII. LOW-ENERGY KINK IN THE 
QUASIPARTICLE SPECTRA 

We can now compare the quasiparticle spectra pre- 
sented in the previous section with the spectra measured. 
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e.g. by photo-emission, in cuprate-based systems. Un- 
like the hnear dispersion predicted by simple band cal- 
culations, in recent years a series of experiments on the 
electronic structure of many H-TC superconductor com- 
pounds has shown sharp breaks in the dispersion of spec- 
tra , " kinks" , at binding energies of the order of 50-80 
mg Y89,90,9i,92,93 ^ Suddcu changes in the quasiparticle ve- 
locity were reported by a factor two or more. This break 
in the dispersion is evident at and away from the d-wave 
node line, and its magnitude shows doping and tempera- 
ture dependence. Kinks may provide useful information 
on the nature of the coupling between electrons and pos- 
sible single-particle or many-body excitations, which are 
at the origin of strongly-correlated many-body properties 
of the system. In H-TC superconductors these feature 
have been associated either with phonon a^^i^° or spin- 
fluctuation based^iiiS^ pairing mechanisms. 

In Fig. [25] the dispersion we observe in the quasiparti- 
cle band of our result around k = (f, f) is similar to the 
experimentally observed kink (see e.g. TeiM^). A neat 
linearly dispersing quasiparticle crosses the Fermi level 
{lu = 0) in the k-uj plot, but around lo — 0.2t, the 
dispersion suddenly changes in slope and gets more inco- 
herent, as evidenced by the broadening spectra. If we set 
t ~ SOOmeV as order of magnitude, we have that the kink 
appears at w 60meV, in good agreement with the ob- 
served experimental energy range. The kink present in a 
wide range of doping, from the under-doped to the over- 
doped regime, and its slope is increasing with decreasing 
doping. 

This new energy scale in H-TC superconductors arises 
the debate on the possible nature of the electron-electron 
coupling. In order for superconductivity to take place 
in metals, it is necessary that electrons bind into pairs, 
which condense in a phase-coherent quantum state. In 
standard BCS superconductivity, coupling between elec- 
trons and phonons (lattice vibrations) drives the forma- 
tion of the pairs. The existence of the kink low-energy 
scale, not explicable in band theory calculation, may give 
an hint on the low-energy nature of the interaction be- 
tween electrons, and hence help revealing the pairing 
mechanism of the unconventional H-TC superconductiv- 
ity. The single-band Hubbard Model studied in this pa- 
per does not take into account phonons by construction. 
Therefore the presence of the kink supports the idea that 
the origin of these features are indeed purely electronic, 
in agreement with the DMFT and CDMFT studies of 
ref i^^'^^ . In Fig. we present the spectral function 
A{k,uj) = —^lmG{k,u}) calculated in our theory, con- 
fronted with experimental ARPES data taken from ref.^"*. 
Fig. [26l shows the spectral function A{k, lu) as a function 
of the energy lu, along two vertical cuts in the first quad- 
rant of the Brillouin Zone in correspondence of the nodal 

and antinodal region (fcyiodjO) ("2'^) ~^ {J^nodi^nod) 
(f,f) and (fca„od,0) ~ (7r,0) (kanod, kanod) ^ (tT, f). 

The system is close to optimal doping 6c ^ 0.08. In the 
left column we display the CDMFT calculation, while the 
right column the experimental data, taken from refi^. 



The nodal quasiparticle clearly show a dispersion, that 
from the Fermi level propagates at negative frequencies 
until LU ^ — 0.2i, where it has a sudden broadening, indi- 
cating a strong incoherence. The antinodal quasiparticle 
shows instead a much flatter dispersion at energy cor- 
responding to the superconductive gap, and the antin- 
odal quasiparticle (which shows to have less weight than 
the nodal quasiparticle) does not lose much coherence. 
This plots have a good resemblance with the experimen- 
tal data on the right column. The comparison between 
the energy scale of our results in the left — t < w < 
and of the experimental data — 300eV< a; < shows also 
that the esteem we have used t ~ 300eV is a reasonable 
order of magnitude. 

The same comparison is also presented in the inten- 
sity color plot at the bottom of Fig. [211 where the nodal 
(left side) and the antinodal (right side) dispersion cal- 
culated with CDMFT are displayed in the top row, while 
experimental plots of ref.^^ are in the bottom row. From 
these diagrams, it is more evident the quasiparticle-like 
dispersion in the nodal point, which compares very well 
to the data from photo-emission. In the photo-emission 
data however at w ~ — 0.075eV the quasiparticle appears 
loosing coherence (effect marked by the loss of red color 
that gets yellow) and at the same time the slope of the 
dispersion changes substantially. In the CDMFT calcula- 
tion instead, as soon as the quasiparticle loses coherence 
(around lu ^ —0.2t), the dispersion stops (i.e. is Lu{k) ~ 
constant), and no spectral weight is present in the region 
at smaller energy (lu < —0.2t). The kink appears in our 
result as a "gap" in the lu vs k spectra. A similar phe- 
nomenon takes place at higher energies, as marked by 
an arrow in the top-left panel of Fig. 1231 where an evi- 
dent gap is present in between the low energy band and 
the lower Hubbard band. This could be possibly asso- 
ciated with the "water- fall" dispersion features recently 
observed in many cuprate materials (see e.g. refi^!^) 
at energies ~ 350 — 600 meV, where the same fc- vector is 
marking the dispersion Lu{k) in a rather extensive range of 
energy. In our calculation this phenomenon is described 
as a gap in the dispersion. The antinodal point, where 
the superconducting gap is present, is non-dispersing at 
low energy (right-bottom side in Fig. [26]), and the re- 
semblance of our result with the photo-emission data is 
rather good. 



VIII. THE HALL RESISTIVITY TO DETECT 
THE TOPOLOGICAL PHASE TRANSITION OF 
THE FERMI SURFACE 

We have shown in section [Vj that in the under-doped 
region (S < 6p ^ 0.06) a different regime sets in our 
solution, marked by a topological transition from a 
large to a small Fermi surface (which reduces to an 
arc or pocket). The strong reduction of the Fermi 
surface area should correspond to a strong reduction 
of the carrier density too. This could be in principle 
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FIG. 26: (Color online). Top panel: The low energy spec- 
tral intensity A{k,u}) vs ui calculated within CDMFT (left) is 
confronted with the experimental result of ref.— in the nodal 
and antinodal regions. Vertical cuts are traced in fc-space 
close to the node and the antinode, as traced on the figures, 
k is expressed in units of (tTjTt). Bottom panel: Color- 
intensity plots of the spectral intensity A{k,uj) in the k — uj 
space in the nodal (left) and antinodal (right) cuts of momen- 
tum space (the same used in the top panel). The color scale 
is set X — 0.30 (see appendix IX)). In the bottom row the ex- 
perimental result of ref."24 are displayed for comparison. We 
refer to appendix 1X1 for the choice of the broadening irj{Lj). 
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FIG. 27: Hall resistivity Rh as a function of doping S. The 
topological phase transition of the Fermi Surface is efficiently 
detected by a rapid change of the behavior in the Hall resistiv- 
ity Kh, which indicates a change of the carrier density. Left: 
our theoretical result (units are in the unitary cell volume Vo 
over electron charge e). Right: the experimental result for 
LSCO extracted from ref.—. 



experimentally detected by measuring the Hall resis- 
tivity Rh ~ — which is directly related to the 
carrier density ric, provided superconductivity can be 
suppressed and the underlying normal liquid extracted. 
This latter is far from being a trivial procedure, which 
has many caveats. A standard method is for example 
to apply a magnetic field, which suppresses super- 
conductivity. Recent publications have addressed this 
problem in BSLCO^ {Bi2Sri,5iLaoA9CuOe+s) 
and LSCO^ {La2-5SrsCuOi), where mag- 
netic field up to 65 Tesla was applied, and in 



Nd 



LS'CO^"(Lai.6-2:Ndo.4Srj;Cu04), where more 



moderate magnetic fields (up to 35 Tesla) can totally 
suppress superconductivity. The results of these ex- 
periences on the Hall resistivity show in fact to be 
compatible with a reduction of the carrier density at 
low temperature in an under-doped sample (as we will 
discuss more in detail in the following). 

In our study we can easily estimate the Hall resistivity 
by using a Boltzmann approach, similarly for example to 
the phenomenological study of ref.^ The Hall resistivity 
Rh can be expressed in terms of the in-plane conductiv- 
ity (in the direction of an applied voltage) a^x and the 
transverse conductivity axy (perpendicular to the applied 
voltage) : 



R 



1 



H 



(39) 



where H is the applied field. We can use the low-energy 
limit (similarly to eq. I30p and extract the normal compo- 
nent of the system (by setting T^ano = in eq. [11]). The 
Green's function is conveniently written on the Matsub- 
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ara axis as: 

, (40) 

where a re-normalized dispersion = [^^ + ReEfc(O)] 
and a inverse temperature re-normalized by the 

quasiparticle residuum Z^, are introduced. We consider 
this liquid as the base on which to start applying a Boltz- 
mann theory. The conductivities can then be expressed 
in the first order in H (see e.g. ref4^): 

<Ta;x = -2e T 2^(Wfe) — tt: (41) 

Here, v'^ = dQk/dki, is v = x,y component of the normal 
velocity, f{x) = l/(e^^ -I- 1) is the Fermi function, r is 
the scattering time, which at T = is, for example, given 
by the impurities, and that, for convenience's sake, in a 
first order approach, we can hypothesize fc-independent. 
In the pure model, like the Hubbard Model, this term 
is absent, so we have to add it in order to simulate the 
finite T = resistivity of a real material. As the Hall 
resistivity Rh is anyway independent of a constant r, 
we do not need to make any special further assumption 
about it. 

The Hall resistivity resulting from our calculation im- 
plemented via a mixed periodization and plug into the 
Boltzmann expression is displayed in the left hand side 
of Fig. [21] as a function of doping 5. The unit used to 
display Rh is the ratio between the unitary cell volume 
Vo and the electron charge e. To give an idea of the 
order of magnitude, taking for LSCO (see e.g. refi^) 
an average lattice spacing Oo ^ 4A = 4 x 10~'^mm 
in the Cu — O2 planes, and oi ^ Soq in the c-axis 
^ ~ 1.2 mm'^/C, which well compares with the ex- 
perimental results extracted from ref.*^ (but similar re- 
sults hold for BSLCO^), and presented on the right 
hand side. In the following, in drawing a parallel be- 
tween our theoretical result and the experimental data, 
one should keep in mind that in our study the "optimal 
doping" (which we can only identify with some degree 
of uncertainty as a maximum in the order parameter 
Sc ^ 0.08 < S < Sp ~ 0.06, see also the discussion in 
the conclusions) is situated at a smaller doping than the 
experimental value 6opt ~ 0.17. 

Rh is positive, as expected by the hole-like Fermi sur- 
face in the hole-doped system. As, in approaching the 
Mott insulator, the localization of particles reduces free 
carriers, Rh is generally expected to monotonically de- 
crease with doping. Two regimes are clearly separable 
in our result (left side of Fig. At (5 - (5p - 0.06 

we observe a discontinuity in the behavior of Rh (which 
slightly decreases instead of increasing) , and for S < 6p a 
sudden increase of Rh, related to the topological phase 
transition of the Fermi surface presented in section |Vl 



marks the reduction in the carriers in approaching the 
Mott state. Such a discontinuity is present and more ev- 
ident in the experimental data (right side of Fig. [27]) . 
The presence of a local minimum at optimal doping 
^opt ~ 0.17 represents a crucial point in the results pre- 
sented in reference^^^ . This behavior is interpreted 
in terms of a quantum critical point, associated with a 
change in the topology from a large hole-like Fermi sur- 
face (realized for 6 > 5 opt) to a small (pocket-like?) Fermi 
surface, realized for 5 < ^opt)- According to their point of 
view, the force driving the superconductivity is related to 
the fluctuations around a critical point. This attractive 
force could overcome the mutual Coulomb repulsion of 
electrons, delocalizing and freeing carriers right in prox- 
imity of the quantum critical point. This fact would orig- 
inate the local minimum in the Hall resistivity Rh (which 
is inversely proportional to the number of free carriers). 
The critical point would correspond then to the opti- 
mal doping (5opt, where the critical temperature Tc is the 
highest. 

Our result shows qualitatively very similar trends, with 
a minimum (even if milder) just at doping 6p, where the 
topological transition of the Fermi surface takes place. 
This strongly supports the analysis of our results and 
the conclusion we have derived, in comparison with the 
experimental evidence on the Hall resistivity. 

In our 2x2 plaquette study at T = 0, however, we are 
not unfortunately able to state if dp is or not the optimal 
doping (which is defined at T = Tc). This point is as a 
matter of facts close to the maximum of the d-wave or- 
der parameter (as evident in Fig. [1]). The study of ref4i 
has pointed out that the exact determination of the op- 
timal doping may depend on temperature. Also, we can- 
not demonstrate that 6p ~ 0.06 is a quantum critical 
point, as we have not clearly identified an order param- 
eter or the divergence, to some order, in the free energy. 
Our study however finds good agreement with the exper- 
imental results of ref] ^°'^^i^^ , and it is not in contradic- 
tion with their conclusions. We can in fact associate the 
discontinuity of Rh with a dramatic re-arrangement of 
the electronic structure in the system, corresponding to 
the topological transformation of the Fermi surface. This 
happening is marked by a fast crossover region (between 
6c ~ 0.08 < 6 < Sp ^ 0.06), where the maximum of order 
parameter is located, the pseudogap appears (at ^ ~ (5c), 
the superconducting gap start decreasing (at S ~ Sp). 
All these facts have striking consequences on spectra and 
transport properties, which we have presented through- 
out the paper. 

More developments beyond the 2x2 plaquette 
CDMFT are needed, in order to be able to reveal if in 
reality Sc = Sp, or rather they represent two different 
"transition" points. It would be then important to clar- 
ify if any of these points are a quantum critical point, 
and their exact connection with the optimal doping Sopt 
of the system. All these are important open questions 
left for the future studies. 
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IX. CONCLUSION 

A. Comparison with the resonating valence bond 
mean field theory 

CDMFT can be viewed as a generalization of the 
earlier slave boson resonating valence bond mean field 
theorj*i22ii2L. It is therefore useful to put our results 
in this context. RVB mean field theories had numerous 
early successes, like the prediction of d-wave supercon- 
ductivity and a pseudogap phase having the same sym- 
metry as the superconducting stated. Both RVB mean 
field theories and our CDMFT study are formulated in 
terms of variables defined in a plaquette. The main differ- 
ences stem from the fact that in CDMFT these variables 
are frequency dependent, and thus they are able to prop- 
erly describe coherence to incoherence crossover in the 
momentum space. In the slave boson RVB mean field 
theory there are two very important notions. The slave 
boson order parameter, which measures "Fermi liquid co- 
herence" , namely the emergence of a well defined quasi- 
particle peak in the spectral function, and the "spinon 
condensation" order parameters, both in the particle- 
particle and particle-hole channel. The position of the 
quasiparticle is shifted by a Lagrange multiplier which, 
adding to the bare chemical potential, forms a quasipar- 
ticle chemical potential. The dynamical mean field pic- 
ture of CDMFT develops these ideas, allowing them to 
acquire non trivial dependence in momentum space. 

To begin with, CDMFT describes naturally a high tem- 
perature state showing poor coherence. As discussed in 
reference^^ and Fig. [7l the peaks in the spectral functions 
in the normal state are broad and incoherent. It is only 
in the superconducting state that we can, for the first 
time, extract the Fermi liquid parameters discussed in 
this work. Coherence appears at low temperature in the 
nodal region more than in the antinodal region, as it is 
clearly seen for example Fig. [71 Furthermore the weight 
of the quasiparticle is different (and larger in most of 
the phase diagram) in the nodal than in the antinodal 
region, as seen in Fig. 1101 At small doping quasiparti- 
cles in the antinodal region cannot be defined in a strict 
Landau-Fermi liquid sense, as a quasiparticle peak can- 
not be identified at the Fermi level anymore. As a matter 
of fact a pseudogap opens in the antinodal quasiparticle 
spectra. In the nodal region instead, quasiparticles are 
always well defined and their weight dramatically reduces 
as the doping is reduced. The slave boson RVB picture 
describes the pseudogap phase in terms of the formation 
of spin-singlets, parameterized by two order parameters, 
one describing correlations in the particle-particle chan- 
nel and the second in the particle-hole channel. A simi- 
lar but more complete description is achieved in terms of 
the dynamical anomalous and normal self-energies. No- 
tice the remarkable fact that the anomalous self-energy 
at low temperatures and at low frequencies has a simi- 
lar order of magnitude as the normal self-energy (Fig. [2] 
and Fig. O, indicating that, as in the RVB mean field 



theory, they might both have a similar singlet-pairing 
origin deriving from the approach to the Mott insulator. 
On the other hand, the normal state self-energy has a 
more complicated angular dependence than the simple 
harmonic (d-wave) fc-dependence of the anomalous self- 
energy. In this sense, the normal and the anomalous gaps 
can also be viewed as competing for the same electrons, 
suggesting that further refinements of CDMFT along the 
lines of ref.^"- are worth being pursued. 



B. The scenario presented by our CDMFT results 

While the dynamical frequency dependence intro- 
duced by DMFT well describes the low-to-high energy 
crossover, as well known from the infinite dimensional 
case^, the momentum dependence introduced by a cluster 
DMFT uncovers a wider spectrum of fc-dependent physi- 
cal phenomena, which turn out fundamental in describing 
the approach to the Mott transition in two dimensions. 

In the first part of the paper (section [TTI we have pre- 
sented raw cluster quantities, which are direct output 
of the CDMFT procedure, but which can be only par- 
tially interpreted in physical terms. Nevertheless we have 
shown that the doping-driven approach to the Mott tran- 
sition takes place via an intermediate regime (rising after 
an "optimal doping region" 5c ~ 0.08 > i5 > ~ 0.06), 
where physical properties depart from the standard pic- 
ture of a BCS superconductor. In particular we were able 
to identify two distinct energy-scales, one associated with 
the anomalous component of the self-energy (see e.g. Fig. 
[5]), the other better enlightened in the local density of 
states (see Fig. [S]), showing different doping dependence. 
In the low doping regime the local spectra also show a 
strong asymmetry in as observed in experiments^^ and 
contrary to the expectations from a BCS superconductor. 

In order to physically interpret the genuine cluster re- 
sults, we have restored the lattice translational invari- 
ance (broken in the CDMFT procedure) by introducing 
a periodization scheme. In section IIIII we have justi- 
fied and compared two possible methods, based on pe- 
riodization of the cluster self-energy and the cluster cu- 
mulant respectively. To this purpose, we have discussed 
the physical properties of our system, taking advantage 
from either the cluster results (supported also by the 
comparison with CDMFT results obtained with QMC 
impurity-solver method s^^'^^i^^ ) and experimental obser- 
vation. We have performed a robust test on our ap- 
proach by reconstructing the local quantities from the 
momentum-dependent Green's function G{k,u}) (see Fig. 
I^O)) . which show to be not far from the correspond- 
ing cluster quantities obtained directly in the impurity 
model. 

By introducing a periodization procedure in 
momentum-space we were able to analyze our re- 
sult in terms of experimentally observable quantities 
in the nodal and antinodal points of momentum space 
(section IIV[) . making contact with recent spectroscopy 
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experiments^!^. In particular, we were able to interpret 
the two energy scale in terms of a pure superconducting 
gap (dominant in the nodes) co-existing with a normal 
component gap (related to the pseudogap of the normal 
state and dominant in the anti-nodal region of momen- 
tum space at low doping). We complete in this way the 
work presented in the short publication of ref4^. 

We have then extended our procedure to describe in 
first approximation physical properties in all momentum 
space (section [V]) . The scenario which results present 
an under-doped state where electronic structure under- 
goes a dramatic re-arrangement. We could associate this 
fact with the appearance of a topological phase transi- 
tion of the Fermi surface (see e.g. Fig. I16p . driven by 
the appearance in momentum space of lines of diverging 
self-energy, fingerprints of the Mott physics. While these 
results may be quantitatively different for different peri- 
odizing methods, in terms of a large enhancement of the 
real part of the self-energies they indicate the same qual- 
itative trends. Enhancement of self-energy Sfe is most 
relevant in the cumulant scheme. The formation of a 
pseudogap in the antinodal region, however, only requires 
a large value of the self-energy and not a strict divergence 
(which in any case can only occur at T = 0). Based on 
this physical idea, at finite doping the antinodal region 
is closer to the Mott insulating phase, while the nodal 
region is closer to the Fermi liquid state (see the band 
structure of Fig. [Ml and the quasiparticle peaks Fig. [5]). 
We achieve in this way a consistent picture of the evo- 
lution of the electronic structure with doping within our 
formalism. 

We were able to explain the strong asymmetry in the 
local density of states, observed e.g. in scanning tunnel- 
ing experiments^. We were able to give a complete de- 
scription of the evolution of electronic spectra with vary- 
ing doping, comparing the doped state with the Mott 
insulating one (see section lyl]) . In particular, we have 
shown that the transfer of spectral weight in approaching 
the Mott transition takes place in a strongly anisotropic 
fashion. We have described the behavior of the low en- 
ergy band close to the Fermi level, which is strongly 
renormalized by the interaction and progressively de- 
stroys by reducing doping, opening first a pseudogap in 
the antinodal region. This latter appears to be present 
already in the parent Mott insulator. We have stressed 
how this phenomenon is connected to the aforementioned 
lines of diverging self-energy. We have described how, 
from this underlying structure of the normal component, 
a d-wave superconducting gap rises. We have also shown 
(in section WU^ how the combination of these effects re- 
sults in spectra which show " kink" features similar to the 
ones observed in the electronic dispersion of many photo- 
emission experiments. Finally we have shown, using a 
simplified Boltzmann approach, how the rising of the 
under-doped regime, associated in our study to a topo- 
logical phase-transition of the Fermi surface, is experi- 
mentally detectable from the doping-dependent behavior 
of transport quantities, like for instance a singularity in 



the Hall resistivity. All these properties are compara- 
ble with experiments on cuprate materials, and the good 
agreement we have found supports our study of the evo- 
lution with doping of the superconducting state of the 
Hubbard Model, offering a self-consistent scenario of the 
approach to the Mott transition in two dimension. This 
is most relevant in connection with cuprate H-TC su- 
perconductors, in which the evolution of the electronic 
structure with doping (and its relation with H-TC su- 
perconductivity) remains a fundamental open question. 
In particular, in recent times, the issue of two (nodal and 
antinodal) gaps^i^ or one pure d-wave gap ^°^'^°^'^°^ has 
been at the center of the experimental and theoretical 
debate^iSi^. In this contest, our CDMFT results on 
the simplest electronic model for H-TC materials show 
a good agreement with the two gaps scenario and, at 
the same time, with many other electronic properties of 
cuprates. 

An important question left open in our study is to de- 
termine the exact nature of the small " transition region" 
between Sp ~ 0.06 and Sc ~ 0.08, where fundamental 
changes take place in the physical properties of the sys- 
tem, in going from the Fermi-liquid of the over-doped 
side into the anomalous liquid of the under-doped side of 
phase diagram. In particular, we have shown that lines of 
diverging self-energy appear at the Fermi level, a pseudo- 
gap opens in the one-particle spectrum and a topological 
phase transition of the Fermi surface occurs (see Fig. [TB)) . 
This small region of doping locates, within the numerical 
precision of our result, the "optimal doping", which we 
identify as a maximum in the d-wave order parameter 
(Fig. [T|) . All these observations point towards indicating 
a tight connection between the physical happening taking 
place in this optimal doping region and the H-TC mecha- 
nism. We are not able to state within our study if behind 
lays a quantum critical point (which could provide the 
binding force for a high critical temperature). This sce- 
nario is actually supported by many theorie o^°^'^°^'^°^'^°^ 
and experimental studies^iii^iiii, and it could be that 
further developments of CDMFT (i.e. increasing cluster 
size) could reveal that the two points Sp and Sc actu- 
ally coincide. Or two distinct points could be actually 
present, and one (or both) could have the characteris- 
tic of a quantum critical point. To reveal this, the right 
divergence in any order of the free energy and an order 
parameter should be clearly identified. 

At present it is not possible to demonstrate that our 
solution will survive in the thermodynamic limit (i.e. in 
the infinite cluster size limit), in the sense that the real 
ground-state of the Hubbard Model in two dimensions 
could be another one of the possible competing instabili- 
ties, like, e.g. stripe ordering or antiferromagnctism. The 
latter, in particular, is expected to be the ground-state 
close to the Mott insulating state. Within a (dynam- 
ical) mean-field approach, however, we can study the 
pure paramagnetic phase, showing that it is a relevant 
phase even if it is not the true ground-state of the sys- 
tem. We leave open the physical question of what terms 
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need eventually to be added in the Hamiltonian in or- 
der to make this state a real ground-state. In order to 
make our first order picture more rigorous, further devel- 
opments are needed. These should involve the cumulants 
directly in the self-consistency condition, as proposed for 
example in refi^ii^^, and should exploit the flexibility of 
CDMFT, which can be formulated in terms of a set of 
adaptive patches in momentum space (on this line see the 
recent work-1-). These extensions, as well as the use of 
more powerful solvers, which could allow going beyond 
a 2x2 cluster, are worth pursuing and are left for future 
studies. 



APPENDIX A: DISPLAYING FUNCTION ON 
THE REAL AXIS WITH ED-CDMFT 

Once the ground-state \gs > of the associate Anderson 
Impurity Model (eq. [8]) has been determined via the 
Lanczos procedure, it is possible to determined the zero- 
temperature Green's function via a second Lanczos step 
(see e.g. the review^). To this purpose, one has to take as 
initial vector cj^\gs > (/i denoting the generic cluster-site 
index), and write the ground-state Green's function in 
a continued-fraction expansion, describing the "particle" 
and "hole" excitations: 



with 

G"=>'<(^„) 



(Al) 



(A2) 



The parameters a" and b" {j = 1... number of Lanczos 
steps, a =>, <) output directly from the second Lanc- 
zos stepii^. In the ED-CDMFT procedure the cluster 
Green's function and the self-consistency equation [3] are 
evaluated on the Matsubara axis. This implies introduc- 
ing a parameter /3 which determines the grid of Matsub- 
ara points cjn = {2n— 1)tt/P, and which plays the role of 
a fictitious temperature (our solution is however at zero 
temperature). We can easily analytically continue the 
Green's functions by replacing in the continued-fraction 
expansion (eq. IA2p 



17] 



(A3) 



where lo is the real axis frequency and rj a small param- 
eter used to display the poles. It is difflcult to know 
a priori the smallest value we can assign to rj. It de- 
pends on the physical problem considered, on the size of 
the associated Anderson Impurity Problem used in the 
Lanczos procedure and on the energy resolution imposed 
in satisfying the self-consistency condition on the Mat- 
subara axis (i.e. the parameter /3). A reasonable guess 
for a lower bound value would be for example 77 « 7r//3. 
Moreover we expect rj to be frequency dependent too. 



The Lanczos better determines the ground state of the 
system, and portrays better low-energy properties. The 
uncertainty in the ground-state vector \gs > propagates 
in the determination of the a" and b" coefficients (which 
is a further Lanczos step), and further propagates in the 
periodization procedures (section IIII[) . The error in the 
second Lanczos step turns out much bigger at higher fre- 
quency, and, if rj is chosen too small, this can create enor- 
mous errors, even breaking casuality (which is instead by 
construction always satisfied in the impurity solver out- 
put, i.e. in the cluster Green's function of eq. [5]). A too 
big value of 77 results however in a poor resolution, which 
may hide important features, above all at small energy. 
In order to be able to display at the same time high and 




FIG. 28: Example of the broadening r;(aj) used to display 
spectral function on the real axis. Left: the low energy range 
— 1 < tj < 1 example of Fig. 1201 Right: the wide energy 
range — 14 < < 14 example of Fig. [T^ For a complete list 
of Tj{uj) see table 1X1 

low energy features, we have introduced for convenience's 
sake in sections IIIII (with the only exception of Fig. [3 
see the corresponding caption), |Vl |Vl] and IVIII a uj de- 
pendent 77, arbitrarily choosing a function-shape which 
quickly separates a low from a high energy range (see 
Fig. Hg): 



(A4) 



The branching between the low and high energy regimes 
takes place in a Fermi-function-like step at uj — loq with 
rapidity l/T^,. In Fig. [28] we report as example the 77(0;) 
adopted to display the local density of states of Fig. [20] 
(on the left), where the low-energy range is — 1 < < 
1 , and (on the right) the 77(0;) of the density plots of 
the Fig. [T9I and [2T| where the full energy scale —14 < 
Lu < 14, including the lower and upper Hubbard bands, is 
considered. In the following we report exactly the values 
of the parameters determining 77(0;) for all the figures 
involved: 
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It is clear that spectral peaks displayed with different 
r](uj) at different frequency present uncomparable heights 
and widths. Our task is not however to make this kind 
of comparison, rather to focalize on the position of the 
spectral peaks and compare heights and widths at differ- 
ent doping (but at the same frequency). This has to be 
kept in mind when analyzing these figures. 







FIG. 29: Color scale code adopted in displaying the color- 
density plots. X is the maximum value of the scale. 

Finally, the color-code we adopted to display the 
density-plots (see e.g. figures [H [HI [21 [21 [Ml 123 is 
shown in Fig. 1291 x is the maximum value of the scale, 
that we have chosen according to the picture (see the 
corresponding caption). If the value of the function dis- 
played is bigger than x, the color remains red. 



APPENDIX B: SUPERCONDUCTING 
BATH-PARAMETERIZATION WITHIN CDMFT 

The general form of the associated cluster-Anderson 
impurity Hamiltonian Tiimp can be written: 

Wimp = ^ Ef^„crcl^C^a + ^ X! '^^'T^Mi + 

+ '^kaal^aka + Vk^'yol^Cf.a + h.c. + 

k(T kfj,a 

+ E K7AA<^ + E <7c,.a,. (Bl) 

Here /i, = l,...,Nc label the sites in the cluster 
and E^^a represents the hopping and the chemical po- 
tential within the cluster, eka is the energy level of 
the (fc, cr) orbital in the free electron bath, Vs rep- 
resent the hybridization hopping amplitude either for 
a particle-desctution/construction and for the singlet- 
destruction/destruction (construction/construction) be- 
tween the impurity-cluster and the bath. In order to 
drive the solutions towards physically interesting regions 
of the bath-parameters space, we have introduced a re- 
duced bath-parametrization^^ which allows to exploit the 
symmetries in the square lattice to gain a better physi- 
cal insight of the Green's function symmetries. Moreover 
using fewer parameters the work required by the min- 
imization procedure is faster and the result simpler to 
interpret: 



n 



imp 



E^ + ^E'^^T^/^i 



E- 

mm' (7 



+a^^a?^ - a'^^a'^^ + a^^a^^ - a^^a^J + h.c. (B2) 

The energy levels in the bath are grouped into multiples 
of the cluster size {N^ = 4) with the labels m = 1, - ■ ■ ,Nc 
and a = 1, 2 such that we have 8 bath energy levels e^^ 
coupled to the cluster via the hybridization matrix V^^^. 
Using lattice symmetries we take V^^^ = V^Smfi and 
= e" . A" represents the amplitude of superconduct- 
ing correlations in the bath. No static mean-field order 
parameter acts directly on the cluster sites^^'^^'*. e", 
and A" are determined by imposing the self-consistency 
condition in eq. [1 using a conjugate gradient minimiza- 
tion algorithm with a distance function that emphasizes 
the lowest frequencies of the Weiss field^ (see eq. [T2]) . 
The reduced form eq. IB2I is in fact a sub-case of the 
more general Hamiltonian eq. IBH and a canonical trans- 
formation connects eq. IB2l to eq. IBll In order to see this, 
it is most convenient to express 7i in a Nambu's form, 
introducing cluster-plaquette spinors: 



= (O'lp ■ ■ ■ ,a4T'"U' ■ 
and recast the Hamiltonian: 



iC4i) 



(B3) 



Hi„,p = *^E* + C/^n^tn^i+ (B4) 

a 

where the Hamiltonian coupling constants are expressed 
my 8X8 matrices (leaving for convenience's sake the 
multi-bath index a implicit): 
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In the more general case t4§^, the bath matrix is diag- 
onal, so that Eg"^^ Skk'SkiT, and Eg'^^^ 0. V^r and 
have generally non-zero elements. In the reduced pa- 
rameterization case instead (eq. [B2j, the bath matrix is 
not diagonal but chosen to mimic a cluster-plaquette, by 
introducing in each multi-bath a the same bath-energy 
Eaa on every bath-site, a next-neighbor hopping and 
a nearest next-neighbor hopping <[,": 
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Moreover superconductive d-wave pairing terms A^ ap- 
pear in the bath: 



E 



while the hybridization between the cluster and the bath 
is only normal = and simplified by connecting each 
site of the cluster /i with only one site fc^ of the multi- 
bath a with the same coupling constant Vaa- 



B2] 

era 



To connect the reduced parameterization eq. IB2I to the 
general eq. IBll it is sufficient to diagonalize the hermitian 
bath-matrix Eg via a unitary transformation^?.: 



and we have 



/ Vac. 
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(B5) 



and S J Sq, = 1 . This is in fact the requirement needed 
to have a canonical transformation which preserves the 
fermionic commutation-relations. In fact if we apply the 
transformation to the vectors: 




^tEUst 



(B6) 



is 

Ek slA - 1 (BIO) 

Finally Vq is subjected to the same transformation: 

vP^SlvP (Bll) 
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